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We study the optical properties of crystals with spatial dispersion and show that the usual Fresnel 
approach becomes invalid near frequencies where the group velocity of the wave packets inside the 
crystal vanishes. Near these special frequencies the reflectivity depends on the atomic structure of 
the crystal provided that disorder and dissipation are very low. This is demonstrated explicitly by a 
detailed study of layered superconductors with identical or two different alternating junctions in the 
frequency range near the Josephson plasma resonance. Accounting for both inductive and charge 
coupling of the intrinsic junctions, we show that multiple modes are excited inside the crystal by the 
incident light, determine their relative amplitude by the microscopic calculation of the additional 
boundary conditions and finally obtain the reflectivity. Spatial dispersion also provides a novel 
method to stop light pulses, which has possible applications for quantum information processing 
and the artificial creation of event horizons in a solid. 
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I. INTRODUCTION 



The problem of optical properties of crystals with spa- 
tial dispersion has remained challenging since the original 
paper of Pekar on the optics of exciton bands . Despite 
considerable effort, the complete theoretical description 
of the optical properties of such systems is still missing 



The nontrivial optical features of crystals with a dis- 
persive dielectric function e(w, k) are based on the fact 
that incident light with a given frequency excites several 
eigenmodes with different wave vectors k. This poses 
the fundamental problem that the Maxwell boundary 
conditions, i.e. the continuity of the electric and mag- 
netic field components parallel to the surface, are insuffi- 
cient to calculate the relative amplitudes of these modes 
and consequently to describe physical quantities, such 
as reflectivity or transmissivity. Since the early work 
of Pekar and Ginzburg § this difficulty was usu- 
ally addressed in a purely phenomenological approach 
by introducing so called additional boundary conditions 
(ABC) for the macroscopic polarization. These ABC are 
motivated physically by the microscopic structure of the 
surface, but the choice of ABC is not universal and may 
be controversial, see Ref. Q and Comments to this pa- 
per. Only the complete solution of the microscopic model 
can determine the dependence of the reflectivity on the 
microstructure unambiguously. 

Such a solution was found recently for the first time 
for the reflectivity near the Josephson plasma resonance 
(JPR) in highly anisotropic layered superconductors d), 
which is an interlayer charge oscillation due to the 
tunneling of Cooper pairs and quasiparticles in highly 
anisotropic layered superconductors ||lO|-[l^. Josephson 
plasma oscillations inside layered superconductor may be 
excited by the light incident to the surface of the crystal 
in the geometries (a) or (b) shown in Fig. 1. The JPR in 



layered superconductors is the simplest example, which 
illuminates the effects of spatial dispersion and the dis- 
crete atomic structure on optical properties in strongly 
anisotropic materials. Here we will describe the method 
of the calculations in |^ in more detail, generalize our 
results for the JPR to different geometries, discuss the 
various transmission and reflection coefficients in a finite 
size sample and point out perspectives to stop light with 
the help of spatial dispersion. We also stress that the 
discrete atomic structure within the unit cell can have 
similar effects as spatial dispersion. 





FIG. 1. The geometry of the layered system showing the 
incident and reflected light at the surface of incidence (a) 
parallel and (b) perpendicular to the layers. Interlayer charge 
oscillations (vertical arrows) are excited by the component of 
the electric field perpendicular to the layers. 

In the framework of the Lawrence-Doniach model 
|l3| (interlayer Josephson coupling) we can describe 
both layered superconductors with identical intrin- 
sic Josephson junctions (such as Tl-2201 |l^ , |l5| , 
Bi2Sr2CaCu208 [E6|, the organic material k-(BEDT- 

TTF)2-Cu(NCS)2ll3'lll (LaSe)(NbSe2) |l|,|§) and 
compounds, where different junctions alternate like in 
SmLai_^Sr^Cu04-5 |2l|-||, Bi-2212/Bi-2201 ^ or 
atomic scale YBCO/PrBCO superlattices Thereby 
we take into account not only the dispersion of the plasma 
mode caused by the inductive interaction of currents par- 
allel to the layers, but also the c-axis dispersion due to 
charge fluctuations on the layers |29|-|3^ 
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The JPR is an ideal choice to illustrate the effect of 
spatial dispersion and the atomic structure on optical 
properties both theoretically and experimentally. First 
of all, recent optical experiments on the layered super- 
conductor SmLai_a;Sra;Cu04_5 with T* crystal structure 
showed evidence that the spatial dispersion of the Joseph- 
son plasmon in the direction perpendicular to the layers 
is important pl|-[25|. For incidence parallel to the lay- 
ers, see Fig. 1(b) at 6* = 0, two peaks at « 7 and « 12 
cm""'^ were observed in reflection, which can be natu- 
rally understood as the JPR |^ of alternating intrinsic 
junctions with SmO or LaO in the barriers between the 
CuOa-layers |3-|2|,||. The very high ratio of the peak 
intensities, about 20, cannot be explained in a disper- 
sionless model Q and it points to a quite strong c-axis 
dispersion of the plasma modes due to charge variations 
Secondly, from the theoretical point of view the 
well established Lawrence-Doniach model |Q formulated 
in terms of finite-difference equations for electromagnetic 
fields and phases of the superconducting order parameter 
is sufficient to provide a complete microscopic description 
and can be solved analytically. Lastly, it is fortunate that 
the damping due to dissipation is low, because at low 
temperatures the JPR frequency is well below the su- 
perconducting gap and the quasiparticles responsible for 
dissipation are frozen out. Otherwise it would strongly 
overshadow the effects of dispersion or the atomic struc- 
ture as described below. 

Extracting the strength of the c-axis dispersion in high 
temperature superconductors is important in its own, as 
the dynamics of Josephson oscillations in layered super- 
conductors is strongly influenced by it ||3^,^,^. It is 
also intimately connected with the electronic compress- 
ibility of the superconducting Cu02-layers, which is hard 
to measure in situ otherwise and contains unique infor- 
mation about the electronic many-body interactions in 
the layers. 

From a more fundamental point of view, we show that 
in the presence of spatial dispersion the conventional 
Fresnel formulas for reflectivity and transmission have 
to be modified substantially near certain frequencies, if 
both the dissipation and the crystal disorder are weak. 
Usually it is assumed that the optical properties of crys- 
tals are completely determined by average, bulk proper- 
ties described by a frequency dependent dielectric func- 
tion e{uj), but not by the explicit spatial dispersion (k- 
dependence) or the specific atomic structure of the crys- 
tal [implicit spatial dispersion). This is based on the 
notion that the wave length of light is much larger than 
the atomic length scales and therefore light is expected 
to be infiuenced only by averaged properties of the crys- 
tal. Here we will stress out that this approach breaks 
down, if the group velocity, Vg — 9a;(k)/(9k, of the wave 
packet of the optical excitation with dispersion ti^(k) be- 
comes small. The physical reason for this breakdown of 
the macroscopic theory is the appearance of a small effec- 



tive wave length, Ag = Vg/u, related to the slow motion 
of the wave packet, which can be comparable with the 
interatomic distance. 

The conditions, when the group velocity becomes 
small, can be most easily seen for an isotropic medium 
described by the dielectric function e{uj,k). Then the 
dispersion relation of an optically excited eigenmode is 
c^fc^ = bj'^e(uj,k) = U!^n'^{u!,k). For a transversal wave 
the implicit derivative of this equation with respect to k 
leads to 

^9 dk kl+^^- 

From Eq. (|l|) it is clear that light can be slowed down 
(a) due to a strong frequency dispersion u!dn{uj)/dw ^ 
1 (as discussed in [ |39| , ^ ), (b) due to a small value 
1 — kdlnn/ dk, i.e. when the spatial dispersion is strong 
or, (c) when the wavevector k becomes large. In the ab- 
sence of spatial dispersion in the dielectric function the 
conditions (a) and (c) are fulfilled at frequencies corre- 
sponding to a pole in e(w), where both dn/duj and the 
wavevector k are large, cf. k'^ oc e(w). Furthermore, it is 
expected that in the same frequency region the dielectric 
function is also quite sensitive to the wave vector, i.e. 
explicit spatial dispersion is significant, cf. case (b). 

Accounting for the wavevector dependence of the di- 
electric function, in general leads to multiple solutions of 
the dispersion relation c^fc^ = u'^e{uj, k) for the wave vec- 
tors kzp, p = 1,2, along the direction z perpendicular to 
the surface at given uj in the geometry shown in Fig. 1(a). 
As it will be derived below, only the light-like modes 
with small \kzp\ contribute significantly to the transmis- 
sion and the usual one mode Fresnel result is recovered, 
if \kzi\ ^ \kz2\- On the other hand, the conventional de- 
scription breaks down, when both \kzp\ are comparable 
and contribute to the optical properties. This happens if 
a pole in the dispersionless theory, which corresponds to 
the cases (a) and (c) of low group velocity, is regularized 
by the introduction of spatial dispersion. 

Depending on the type of the spatial dispersion the 
excited modes may be both real (propagating modes) 
or one wave vector may be real, while the other one is 
complex (decaying mode). This leads to two types of 
critical frequencies, where the Fresnel approach becomes 
invalid. Namely it occurs at frequencies uje, where both 
kz are real and \kzi \ ~ \kz2\, and at frequencies uji, where 
kzi ~ ikz2- 
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FIG. 2. Schematic mixing of a transverse optical phonon 
characterized by anomalous dispersion with a propagating 
electromagnetic wave leads to an extremal point uje in the 
lower polariton band, where the group velocity vanishes. Just 
below the frequency uj^ two modes with similar wavevectors 
propagate. 

When both modes are propagating, Vg vanishes at fre- 
quencies LOf. due to strong spatial dispersion, the case (b) 
mentioned after Eq. (Q), see Fig. ||. In general, this case 
occurs, if the eigenmodes of the crystal, when decoupled 
from electromagnetic waves, have a dispersion opposite 
to that of the electromagnetic wave. Generic examples 
are phonon modes with anomalous (decreasing) disper- 
sion mixing with propagating light of normal dispersion, 
which form a polariton (cf. Fig. ^), or the Josephson plas- 
mon with normal dispersion interacting with screened 
electromagnetic waves in a superconductor, which show 
an anomalous dispersion, see Sects. III.B and IV be- 
low. As the main consequence, near frequencies uje the 
transmission coefficient into the crystal is not determined 
solely by the dielectric function, but crucially depends on 
the microstructure of the crystal near the surface, if both 
dissipation and disorder are very low and the system is 
strongly anisotropic. We will also show that interfering 
multiple propagating waves create a behavior similar to 
intrinsic birefringence and affect strongly the transmis- 
sion through the crystal and multiple reflection. 

In the second situation (one mode is propagating, while 
another is decaying) the Fresnel approach breaks down 
near frequencies w^, where the moduli of the wave vec- 
tors of two excited modes become equal. Near these fre- 
quencies both \kzp\ become large, of the order of the in- 
verse interatomic spacing, which leads to a small, but 
finite group velocity Vg as described in the case (c) after 
Eq. (|l|). This occurs, for example, for Josephson plas- 
mons with anomalous dispersion in a crystal with differ- 
ent alternating junctions, where one plasmon has normal, 
while the other one has anomalous dispersion, see Sect. 
IV below. As near the frequencies uji only a single mode 
propagates into the crystal, the transmission coefficient 
is significantly suppressed in comparison with resonances 
at extremal points We, where incident light excites two 
propagating modes. 
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FIG. 3. A pole in kliuj) in the case without spatial 
dispersion (dashed line) indicates the importance of small 
length scales due to the low group velocity (cf. Eq. (jl|)) 
and the breakdown of the macroscopic theory based on a 
fc^-independent dielectric tensor. For an isotropic system 
this corresponds to a singularity in the dielectric function, 
kl oc e(cij) This pole is regularized, when spatial dispersion is 
taken into account, and depending on the sign of dkl/du an 
extremal point We, where the group velocity vanishes, appears 
as shown in a), or as shown in b), the singularity transforms 
into a special frequency tOi, where fc^i — 



'kz2- 



In Fig. ^ it is demonstrated schematically, how the 
critical frequencies uje and w^, where the amplitudes of 
the exited multiple modes equal, \kzi\ = \kz2\, develop 
from a singularity in the one mode theory, which neglects 
the fcz-dependence of the eigenmodes. In the simplest 
case of an isotropic medium, which was considered after 
Eq. (|l|) , the dispersionless dielectric function and squared 
wavcvector amplitudes are proportional, e(a;) oc fc^, and 
their poles coincide. The breakdown of the one mode 
Fresnel theory at these points is already anticipated from 
the low group velocity Vg, due to the large frequency 
dispersion, \dkz/du}\ S> 1, and the large \kzp\ near the 
pole, cf. case (a) and (c) in the discussion after Eq. (^. 

If for the crystal dispersion (without coupling to elec- 
tromagnetic waves) dk^/ dio < 0, an extremal point uJe ap- 
pears below the singularity and at this frequency lj,, the 
group velocity Vgz — doj / dkz vanishes and two propagat- 
ing modes with kzi — —kz2 are excited, see Fig. I (a). In 
a similar way, at the extremum of uj{kz) above the singu- 
larity the imaginary excited modes merge, k'^^ = k'^2 < 0> 
while in the intermediate frequency region the solutions 
k'^p are complex. On the other hand, if dkl/duj > 0, 
the singularity in the dispersionless one mode theory is 
transformed to a special point cui, where the amplitude 
of the excitation equals, but one is propgating and the 
other decaying, k"^^ — — fc^j- 

Remarkably, a special point uji can appear, when the 
group velocity is small, even without a wavevector de- 
pendence (i.e. without explicit spatial dispersion) in the 
dielectric function due to the atomic structure in the 
unit cell alone {implicit spatial dispersion). Generally for 
each polariton band a real or imaginary mode is excited, 
but usually inside one band the second wave associated 
with the off-resonant excitation of the other bands can 
be neglected. Here it will be shown that this assumption 
breaks down, when the group velocity becomes small. 
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e.g. for large amplitudes of the wavevectors, cf. case (c). 
Thereby the system with alternating plasma resonances 
like SmLai_a;Sr3;Cu04_5 with light incident parallel to 
the layers (Fig. |l|(b) at 9 = 0) presents a generic ex- 
ample, as in this case the wavevector kz perpendicular 
to the layers (explicit spatial dispersion) vanishes due to 
the homogeneity of the incident beam. In a macroscopic 
theory the electrodynamic response to the electric field, 
which is averaged within the unit cell, is decribed by the 
effective (average) dielectric function eduj), 
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ec(w) 



1 



1 



eci(w) ec2(t^) 



(2) 



Thereby a pole in ec{uj) appears between the zeros of 
eci{uj) = eco(l - Lolo i/uj^) {I = 1,2, Cco background di- 
electric constant), which correspond to the plasma fre- 
quencies iOcoj in the different junctions This in- 
dicates the breakdown of the one-mode Fresnel approach 
and the necessity to account properly for the second solu- 
tion. Obviously, similar consequenses of such a "discrete" 
implicit spatial dispersion are expected generally for any 
crystals with multiple optically active crystal bands of 
the same symmetry. 

Both the behavior near and tOi are in contrast to the 
conventional Fresnel theory and to the common believe 
that spatial dispersion of crystal modes or the atomic 
structure do not create measurable effects of order unity 
in optical properties, but only enter in negligible correc- 
tions proportional to the ratio of atomic scales and the 
wavelength of light. In fact, the Fresnel results have to 
be modified significantly in a narrow interval near the 
frequencies and coi, but only in perfect crystals with 
very weak dissipation. 

Finally, we point out that the vanishing of the group 
velocity at extremal frequencies oje due to spatial disper- 
sion of the crystal modes provides a novel way to stop 
light pulses dynamically. Recently it attracted a con- 
siderable interest to diminish the light velocity strongly 
with the help of frequency dispersive gaseous media as de- 
scribed by case (a) after Eq. (0). From a practical point 
of view, our suggestion based on the k-dependence of the 
dielectric tensor allows to use slow light in a solid state 
device for the processing of information. In particular, 
the sensitivity of the group velocity in solids to the ex- 
ternal fields could be used to store quantum information 
in the form of photonic qubits, as required for optical 
quantum computers ]4l| ]. Our novel solid state proposal 
to stop light might be of advantage compared with the 
realizations using gaseous media, as it is easier to scale 
to larger system sizes and more complex devices. By ad- 
justing an inhomogeneous external parameter, like the 
magnetic field for the JPR, a spatially inhomogeneous 
profile for the group velocity can be imprinted. Such 
conditions can simulate in the laboratory the behavior of 
light in a curved spacetime, as realized in astrophysical 



situations, e.g. near the event horizon of a black hole 

Previously the spatial dispersion of the Josephson 
plasma mode and its effect on the propagating electro- 
magnetic waves in layered superconductors with identical 
Josephson junctions was discussed by Tachiki, Koyama 
and Takahashi |^l[|. They realized that the mixing of 
plasma modes with electromagnetic waves can lead to 
two propagating waves with different wave vectors for 
the same frequency. However, implications of this fact 
to the optical properties, like reflectivity, were not dis- 
cussed. Van der Marel and Tsvetkov present an ef- 
fective dielectric function for the system with alternating 
Josephson junctions and charge coupling within the unit 
cell for the special case of incidence parallel to the layers, 
but they did not account correctly for the dissipation due 
to the conductivities and for the nontrivial effects of the 
"discrete" spatial dispersion mentioned above. 

The paper is organized as follows: In the first part, 
we derive in general the optical properties of an uni- 
axial crystal with explicit spatial dispersion along the 
symmetry axis in the dielectric function using additional 
boundary conditions with one phenomenological param- 
eter (Sect. II). In the second part, we confirm these re- 
sults for oblique incidence in the microscopic model for 
the JPR accounting for the atomic (layered) structure. 
Thereby the ABC are derived and analytical solutions 
for systems with identical (Sect. Ill) and two different 
alternating (Sect. IV) Josephson junctions are obtained. 
In Sect. ^ the atomic structure is taken into account 
to derive the reflectivity in the incidence parallel to the 
layers. Technical details are given in the Appendices. 



II. MACROSCOPIC APPROACH FOR 
CRYSTALS WITH SPATIAL DISPERSION 

In this Section we derive the dispersion relation from a 
macroscopic dielectric tensor (Sect. II. A), calculate the 
transmission coefficients into (II. B) and through (II. C) 
the crystal using a phenomenological ABC and close with 
some further remarks, concerning e.g. future applica- 
tions, like the stopping of light (II. D). 



A. Dispersion relation 

We consider the geometry of the incident and reflected 
light as shown in Fig. |^. The wave vector of the in- 
cident light with frequency io for the geometry shown 
in Fig. |l|(a) is ko = {uj sm0/c,O,Lu cos0/c), while for 
Fig. |l|(b) ko = {Locos6/c,0,ujsm9/c), where the z-axis 
is perpendicular to the layers (it coincides with the c- 
axis of the crystal). The incident (quasi- monochromatic) 
electromagnetic wave is assumed to be P-polarized, i.e. 
the electric field E(r, t) = E{uj, k) exp(ikr — iwt) is in the 
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plane defined by ko and the normal of the surface {xz- 
plane) , while the magnetic field B has only a component 
in y-direction. S-polarization is not considered here, as 
an electric field parallel to the layers does not excite the 
JPR studied below. 

In the macroscopic approach used here we describe the 
crystal by a dielectric tensor, which is averaged on atomic 
scales within the unit cell, but can depend on the wave 
vector (explicit spatial dispersion), and study the effects 
of the intrinsic microstructurc (implicit spatial disper- 
sion) in Sect. V. 

In the following we will consider highly anisotropic uni- 
axial (layered) crystals with the dielectric function com- 
ponents ec(w, /Cz) along the c-axis (z-axis) and ea(<^) in 
the ab (xy) plane along the layers in a parameter regime 
appropriate for the JPR. In edi^Tk^) we account for a 
collective mode (JPR in our case), which is strictly lon- 
gitudinal with the dispersion uJc{kz) for = 0, i.e. 
ec[uj = ujc{kx — 0,kz),kz] — 0, and whose polarization 
is mainly in the c-direction for any k^ due to the strong 
anisotropy, \ea\ ^ led, near the JPR. We neglect the 
eigenmode, which is polarized parallel to the layers for 
fca, = 0, as it is of much higher frequency than the JPR. 

From the bulk Maxwell equations for the Fourier com- 
ponents, 



ck^By = -iveduj, k^)E^, 
k^E^ - k^Ex = -{uj/c)By, 
ckzBy = ujeaiuj)Ex 

follows directly the dispersion relation. 



kl 



ec(w,fc^) ea(w) c2 



(3) 
(4) 
(5) 



(6) 



of the eigenmodes in the crystal. 

For the geometry shown in Fig. |l|(b) and neglecting 
the discrete layered structure in z-direction, we obtain 
analogously k^ — koz = ujsinO/c of the excited crystal 
mode, while the dispersion relation, Eq. (||), gives a sin- 
gle solution for fc^. Hence, the usual Fresnel description 
is generally valid, except where l/cj-l becomes large, e.g. 
at the poles of ec(w), see Eq. (||). At these points the 
implicit spatial dispersion due to the atomic structure in 
the unit cell in multiband systems has to be taken into 
account. Then multiple solutions fc^ of the dispersion 
relation contribute, which will be discussed for the JPR 
with alternating junctions in Sect. V below. 

In the geometry shown in Fig. |l|(a) we obtain from 
the translational invariance parallel to the surface the 
wave vector component k^ = \iax = tosinO/c, and the 
dispersion relation determines the solution(s) for the z- 
component kz{Lu,9) of the modes excited by the incident 
wave. 

In a crystal described by the dielectric functions 
ea,c(w), which is independent of the wave vector k, the 



dispersion relation Eq. (^ has a unique solution k1{uj). 
The Maxwell boundary conditions (MBC), requiring the 
continuity of the parallel components Ex{z) and Hy{z) at 
the surface z = 0, immediately give the Fresnel formula 
for the reflection coefficient R = |rp and the transmis- 
sivity T = 1 — R into the crystal. Here 
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Exjz ^ 0) 
By {z = 0) cos ( 



(7) 



When in a highly anisotropic crystal the eigenmode 
with electric field approximately parallel to the layers is 
neglected, the effective dielectric function is given by 



'foff 



no 



ea{uj) cos 9' 
where the refraction index is 

no = ckyX^)/bJ = y'eo(w)(l - sin^ 6'/ec(w)). 



(8) 



(9) 



This suggests that for an anisotropic crystal in this geom- 
etry the critical frequencies, where the refraction index 
no becomes large and the Fresnel theory breaks down, 
appear at zeros of tc{^) rather than at poles of the di- 
electric function, as for isotropic system discussed in the 
introduction (cf. Eq. (^) and Fig. ||(b). 

If the dielectric function, ec(w, kz), is dispersive in the 
c-direction, Eq. (^) has multiple solutions for k1{uj) [§,0. 
In the following we restrict ourselves to the simplest case 
of four (in general complex) solutions ±ni,±n2 for the 
refraction indices. 

Generally, in a crystal of finite thickness, where the 
(multiple) back reflection from the second surface is taken 
into account, all four solutions ±ni^2 have to be consid- 
ered. For simplicity, we will consider in the following 
mainly a semi-infinite crystal in the half space z > 0, 
where only two of the solutions are physical. When dis- 
sipation is low, for quasi-monochromatic wave packets 
the direction of the energy transfer is determined by the 
Poynting vector S, which is oriented along the group ve- 
locity Vg =aw/ak [|: 



(10) 



w = 



IQtt 



duj 



e„ez 



duj 



■e,e: 



Here W is the high frequency average of the energy den- 
sity. In agreement with the causality principle the group 
velocity of propagating modes in the c-direction, Vgz = 
duj{kz)/dkz, should therefore be positive. Note that in 
the case of normal (anomalous) dispersion this requires 
the real part of the wave vector kzp (modes p = 1,2) and 
of the refraction index np = ckzp/uj to be positive (nega- 
tive). When dissipation is taken into account, this rule is 
equivalent to the condition that the eigenmodes should 
decay inside the crystal, i.e. lm.{kzp) > 0. 
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This has in particular consequences at extremal fre- 
quencies Ue of the dispersion relation Re(a;(fcz)), where 
the group velocity Vgz = vanishes and two branches, one 
with normal and another one with anomalous dispersion 
merge, see Fig. ^ At these points the two solutions for 
kz , which are real in the absence of dissipation, have the 
same amplitude \kz\, but different signs, 

Re[ni{uje) + n2{uj,)] = 0. (11) 



B. Transmissivity T on surface 

In the macroscopic approach the electric field Ez and 
the polarization Pz in a semi-infinite crystal with a single 
atom in the unit cell and with the background dielectric 
constant eco can be expressed as 



Ez{z) ^ X! Ezikzp)e:iq){ikzpz), 



(12) 



0=1,2 



Pziz) = ^ Ez{kzp)xc{kzp)exp{ikzpz), (13) 

p=l,2 

47rxc(fc^) = ec(w, kz) - Cco, (14) 

while the equations for E^ and By, which enter in Eq. (|^), 
are similar. In order to determine the amplitudes Ez (kzp) 
of the different eigenmodes we use the most general ABC 
proposed by Ginzburg |S| 



Pz{z) + i{dPz/dz) = 0, z ^ 0, 



(15) 



where the length scale £ is a, phenomenological param- 
eter to be determined from the microscopic model. In 
systems with inversion symmetry we can use Xc(w, fcz) — 
Xc(w,0) kz for kz ^ and obtain 



J2Ezikzp){l+i^n.p)^0, ^^Lue/c. 



(16) 



in leading order in ta/n^ 1 and l/|nin2| <C 1. In this 
limit Eq. ( p^ ) and the following results are confirmed 
microscopically for the JPR in Sect. Ill and IV, while 
in general corrections involving field components parallel 
to the surface have to be considered in Eq. (p^. Using 
Eqs. (^), (^ and (p^), we derive (near the resonance) 
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711^2 



(17) 



We see that in the case of multiple eigenmodes in the 
crystal the optical properties like the refiectivity gen- 
erally cannot be expressed by the refraction indices Up 
alone, which are determined by the bulk dielectric func- 
tions ea,c via Eq. (^, but also depend explicitly on the 
parameter ^ introduced by the boundary conditions. 

As the wavelength A of light is larger than all length 
scales related to the atomic structure of the crystal or 



to the change of the polarization at the surface, we can 
assume ^ ~ £/A <C 1. Therefore the term ^nin2 can be 
neglected everywhere except at the extremal frequencies 
We, where Re(ni -I- ^2) — 0. 

If in addition the amplitude of one excited mode is 
large, i.e. |rt2| ^ |"i| and |ni?i2| ^ \ea\, the conven- 
tional one mode Fresnel result, Eq. (||), is obtained for 
the mode with smallest n. In Fig. |^ it can be seen that for 
the phonon polariton away from the extremal frequency 
LUe this condition is fulfilled and only the usual light-like 
mode remains. 

Deviations from the usual Fresnel theory are therefore 
expected, when the amplitudes of ni and 712 are compa- 
rable and both modes play a role. The resonances in the 
transmissivity are located in these two mode frequency 
regions and we distinguish the cases that (i) both ex- 
cited modes are propagating (ni^2 real) or (ii) one mode 
is propagating, while the second is decaying (m real, 712 
imaginary) . The appearance such type of special frequen- 
cies We, where ni = —712, and Wj, where 7ii = in2, near 
a pole in the refraction index of the dispersionless one 
mode theory is schematically shown in Fig. ^ (the index 
of uji reminds of the factor ±z between the solutions 771^2)- 

(i) For two real modes 771^2 we have Re(77i + 712) = 
at the extremal point lu = oje, when causality is taken 
into account, see Eq. (|Tl|). Then, if the dissipation is 
weak in addition, e.g. Im(77i + 772) ^ |C"-i'^2|7 only the 
term i^nin2 in Eq. ( p7| ) remains, ^(we) is imaginary and 
T(tt'e) = 0. The transmissivity T reaches its maximum at 
the frequency Wc,max slightly above iVe- At this frequency 



(ni + 712) = Ea ^^^-l^^2(cos 

re,„,ax = 2/[(l+e'eaCOs2 



+ 



(18) 
(19) 



It is pointed out that both the position Wc.max of the 
resonance in i? or T and its amplitude are determined 
not solely by the imaginary part of Sa.c as in the disper- 
sionless case, but also by the surface parameter ^. This 
correction is important for highly anisotropic systems, 
where ^e^ ^ 1, although ^ <C 1, as it is realized for the 
JPR (see Eq. (^)). We see that in the absence of dissi- 
pation To^max depends on ^ and is generally smaller than 
the Fresnel result Tmax = 1, see Fig. |[ Physically this 
result refiects the fact that the low group velocity near uje 
introduces a small length scale Ag = ^g/i-^i which makes 
the variation of the polarization Pz near the surface rel- 
evant and indicates the breakdown of the translational 
invariance on the atomic scale I. Note that the oppo- 
site signs of the refraction indices ni_2 near We due to 
causality are essential for the dependence of Tcmax on ^. 
The vanishing of ni +n2 at lUq (see Eq . (pl| )) in Eq. (ITT 
and its consequences in Eqs. 



noted previously 



and (M) have not been 

We also note that the results 

in Eqs. ( p_8D and (If ) cannot be obtained from the ABC 
proposed by Pekar | yj^] , which neglects the derivative in 
Eq. (|l|). 
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(ii) In the case, when ni is real, while 77.2 is imaginary 
without dissipation, we anticipate that T is strongly sup- 
pressed, because both modes are excited by the incident 
light, but only a single mode propagates into the crys- 
tal. This situation occurs e.g. in superconductors when 
the dispersion of the collective mode is anomalous (cf. 
Fig. |l| in Section |^). T{uj) is in this case peaked at 
critical frequencies LOi near uJc^ where n2 = —ini with 
ni < 0. Here we obtain for the maximal transmission 
coefficient 



Ti, 



2m 



e„ cos ( 



(20) 



so that T{u)i) <^ J^^inax- This difference in the reso- 
nance amplitude, depending whether two or one prop- 
agating modes are excited, cannot be described in the 
one mode Fresnel approach without spatial dispersion, 
where in both cases a single propagating mode is excited 
and the transmission amplitudes are comparable. This 
observation and the strong deviation from the conven- 
tional Fresnel result is confirmed below for the JPR in 
Sect. IV. In contrast to the situation (i) near extremal 
points LJe, the parameter ^ is irrelevant near w^. 



C. Transmission through thin film 











\ Ri/ / P.. 


t, ' 


lUi 7 / 
















FIG. 4. Transmission amplitudes of the wave with refrac- 
tion index rip into (tp) and out of {jp) the crystal and multiple 
reflections at the first [Rn) and the second surface {Ppp')- 

We now study the transmission and back reflection 
of the multiple excited modes in a thin film of finite 
thickness L, see Fig. ^. For the ratio of the magnetic 
field tpB™ of a partial wave with the refraction index Up 
[p — 1,2) excited in the crystal to that of the incident 
wave B™ we obtain 



2(1 - i^np) 
C("2 - ni){\ -f k) 



(21) 



We will see that \tp\ > 1 for JPR, e.g. the fields of the two 
partial waves are enhanced, but have opposite direction. 
Note that the transmissivity T follows from the ratios of 



the z-components of the Poynting vectors, Eq. (pT|), and 
that T ^ \ti -ht2p- 



In2l»lnil 




FIG. 5. Reflectivity R — 1 — T and ratio \Tp\ of the outgo- 
ing magnetic fields at the second surface of the crystal near 
an extremal point uie with (solid line) and without (dashed) 
spatial dispersion without dissipation (schematically). Com- 
pared with the conventional Fresnel formulas the plasma edge 
in R is at the higher frequency tOe and the amplitude of the 
resonance at aje.max is damped due to ^ in Eq. (^^. The am- 
plitude ~ |Tp(a;e,max)| of the outgoing waves (cf. Eq. (|2^)) is 
strongly suppressed in the frequency region, where R is mini- 
mal, e .g. where the transmission T into the crystal is maximal 
(Eq. (§)). 

At the second surface of a crystal the arriving wave 
with index Tip (p = 1, 2) and the magnetic field amplitude 
By^p creates a wave, which is emitted out of the crystal. 
Its wave vector is kg and we denote its magnetic field 
by TpBy p. Each wave Up also excites two waves with 
refraction indices Upi and magnetic fields Ppp/By^p, which 
are reflected back into the crystal. The ABC Eq. (|l^) at 
z = L for these three waves give 

(1 + i^np)E,,p + (1 - i^m)El% + (1 - iin2)El% = 0, 

(22) 

where E^^p and El^^^, are electric field components at 
the second surface at 2: = L of the arriving and back- 
reflected waves, respectively. We find in leading order in 
C (p = 3-p) 



2np{ni + 712) 



(rii + n2)ea cos 9 + Ea + nin2 ' 



(23) 



- / 1 ^P ("1 + "2) eg cos6>(n2 - rti) + eg ~ nin2 , , 



Ppp 



{n2 — ni) ea cos 9{n2 -f ni) + Ca + nin^ 

{-iy2np{nl - Cg) 

(n2 - ni){ea cos9{ni + 712) + eg + nin2) 



(25) 



At the frequency Wc.max, where the transmissivity T into 
the crystal is maximal, the transmission. 



'^p (^c,inax) 



271,- 



e„ cos ( 



(26) 
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is strongly suppressed in comparison with the conven- 
tional Fresnel result (cf. Fig. At the same point We.max 
the back scattering takes place almost completely into 
the same eigenmode, pn ~ — P22 ~ — 1 + 0{n/ea) and 
P12 — P21 ~ l/[£aC cos 61] ^ 1, while at oj — uje we obtain 
\Ppp\ ^ P12P21 in the presence of spatial dispersion. 

The two eigenmodes of the same polarization interfere 
inside the crystal and for the total transmission Ttot — 
l^totP through the sample we obtain near an extremal 
frequency uje 




ttot = ^ Un e^p{iniUjL/c) 
I 

- [1 + (1 - 2v) cos(2nwL/c)]/2 , 



(27) 
(28) 



where v « (n^/2)^. Therefore, the transmission coeffi- 
cient has oscillatory behavior as a function of the fre- 
quency u! and the sample thickness L due to the interfer- 
ence effect, even if the back reflection into the sample is 
irrelevant. Near the frequency CLie.max multiple reflection 
leads to 



Tt, 



\tt. 



1 



(1 — 2v) cos{2nujL/c) 



1 + p2 _ 2pcos{2nujL/c) 



(29) 



with p = P12P21- 

The difference with conventional birefringence lies in 
the fact that all waves have the same P-polarization. This 
type of so-called intrinsic birefringence has also been ob- 
served in semiconductors for certain directions of prop- 
agation (cf. [Q and references therein), while in the 
present case it appears for an arbitrary angle of (oblique) 
incidence. Alternatively, the effect of spatial dispersion 
can be observed by the splitting of a spatially focused 
incoming beam into two outgoing ones, corresponding to 
the two different group velocities in the crystal (angle 
between rays ^ 10^^ degrees for JPR). 



D. General remarks 

Some additional remarks to the macroscopic approach 
are in place: 

(1) It is pointed out that even if the last term ^ ^ 
in the denominator in Eq. (17) can be neglected for fre- 
quencies far from the band edge near Wc,inax or due to 
dominant dissipation, the interplay of the two modes 
with indices ni^2 can lead to unconventional effects, like 
intrinsic birefringence (Eq. (p8|)) or the suppression of 
the transmission near LUi in comparison with the Fres- 
nel result (Eq. (pO|)). Only in the Hmit \n2\ ^ \ni\ and 
|'^i'^2| S> \ea\ the smallest refraction index determines k, 
Tp and Pppi and the usual one-mode Fresnel description 



is recovered. 

(2) Thereby the existence of a pole in the effective di- 
electric function in the one-mode Fresnel approach is an 
indication of the existence of a special point uJe or uji, see 



and the microcopic confirmation in the Sects. IV 
However, we point out that without further inves- 
tigation of the spatial dispersion or the atomic structure 
these two cases cannot be distinguished. The guiding pic- 
ture in Fig. ^ and the microscopic results for the JPR 
in oblique incidence in Sect. Ill and IV and for phonon 
polaritons Q suggest that special points of type uJe (tOi) 
appear, if light is mixed with a crystal mode of oppo- 
site (same) dispersion. This is seen in Fig. ^ where the 
mixing of the plasma band in the lower (upper) band 
with normal (anomalous) dispersion with decaying light 
creates a special point of type uJe i'-^i)- In Sect. V it is 
shown that special frequencies uii, where = —n^, can 
appear near the pole of the effective dielectric function 
even without k-dependence due to the discrete atomic 
structure within the unit cell. 

(3) It is stressed that the Kramers-Kronig relations ex- 
pressing causality (and sum rules following from them) 
are still valid in the two-mode regime for physical re- 
sponse functions like the reflectivity R or for the effec- 
tive dielectric function tcff — i^"^ extracted from i?, but 
do not apply to the refraction indices Up of the partial 
waves independently [p p5|j4^ . 

(4) We note that beyond the universal electrodynamic 
effects studied above there might also be the necessity 
that the ABC reflect the change of the internal structure 
of the crystal excitations near the surface. This problem 
has been studied in detail for the Frenkel exciton, which 
is quite extended on the atomic scale and whose wave 
function is consequently modified near the surface, see 
Rcf. and references therein. Due to the focus 
on the microscopic derivation of the exciton modes and 
despite a considerable effort, some of the crucial general 
features discussed here have been missed for that system, 
namely the importance of the atomic structure (param- 
eter ^) and the correct causal choice of the eigenmodes 
in a semi-infinite crystal near the extremal points, e.g. 
m -I- ^2 ~ for w « UJe, see Refs. 

In the case of the JPR the effect of the surface on 
the internal structure of excitations turns out to be very 
weak, because the excitations are confined between layers 
on the atomic scale and in highly anisotropic layered su- 
perconductors the layers near the surface are practically 
the same as those inside the crystal. Therefore and be- 
cause we discuss this system only as a generic example for 
general electrodynamic features, which are relevant for a 
large class of systems, we will not address this question 
in the following and assume a dielectric response function 

e,(z,z') = e(z)e(z')e(z-z')- 

(5) The dispersion and the group velocity of phonon 
polaritons has been measured directly by exciting locally 
a wave packet and detecting the time of propagation to 
a separated probe position in the crystal |^^. Future 
experiments of this type with high resolution for long 
wavelengths could also show the existence of extremal 
frequencies We, where the group velocity Vgz vanishes at 
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a finite wave vector as shown in Fig. |[ 

(6) We now comment on the perspectives to stop light 
using spatial dispersion at extremal frequencies uje (cf. 
Fig. H) and compare this method with the alternative 
one, which uses the frequency dispersion of the dielectric 
function |0|. 

The effect of the frequency and/or spatial dispersion on 
the group velocity has already been discussed as a guiding 
principle for an isotropic medium, see Eq. (|l|). In the 
scattering problem depicted in Fig. |l| (b) the component 
kx of the wavevector and the group velocity parallel to 
the layers is fixed by the boundary condition. The signal 
velocity Vgz in z-direction in the anisotropic case [ua = 
^ = follows from Eq. (||), 



duj 
dkz 



d In iig 



a In 



91nne 2 h.2 ' 



(30) 



In the phenomenon of electromagnetically induced 
transparency (EIT) , which has recently been used to cre- 
ate ultra-slow light atomic levels are pumped op- 
tically in such a way that the medium exhibits a sharp 
absorption line in Im(e(u;)) near a resonance frequency wq 
for propagating light. According to the Kramers-Kronig 
relation the frequency dispersion dn/duj of the real part 
of e is therefore quite large, which suppresses the group 
velocity in Eq. (|l|). Spatial dispersion cj(fc) is discussed 
here for the first time as a tool to stop light, although 
a finite drift velocity of a (gaseous) medium has been 
interpreted in this way p^ . 

This effect might be used to realize certain phenomena 
connected with ultra-slow light in a solid, such as the op- 
tical Aharonov-Bohm effect in rotating media or the 
enhanced two-photon interaction via a phonon mode [ pO| , 
which has possible applications in quantum information 
processing. 

Apart from this, the variation of the band structure 
and thus ea,c(r) on scales, which are large compared with 
the wavelength A of light, allows to manipulate the geo- 
metrical optics of light in a solid in a rather simple way, 
e.g. via a space dependent external magnetic field for the 
JPR or pressure for phonon modes. Similar feature have 
been used recently for creating artificially local space- 
time geometries, which are reminiscent of cosmological 
phenomena, such as black holes: e.g. in supcrfluid '^He 
pl| , inhomogenously pumped media with EIT , flow- 
ing dielectrics or solids jssj. In particular, it is possi- 
ble to create a space dependent group velocity profile for 
a given frequency, where Vgz vanishes on some manifold 
in space. At this point the behavior of light is expected 
to be similar to the one near an event horizon of a black 
hole, see 1 42 . Thereby the description in terms of a di- 



electric function breaks down at short length scales. 

From an application point of view, the modification of 
the band structure with the help of an external param- 
eter, opens the perspective to store light pulses dynami- 



cally. Thereby in an ideal crystal the phase information 
of the light pulse or the single photon is stored coherently, 
which makes the device potentially useful in quantum in- 
formation processing ]4l| ]. The limiting factor is clearly 
the decoherence due to disorder or dissipation induced 
by a finite conductivity. For the JPR in Bi2Sr2CaCu208 
the intrinsic decay time due to ohmic losses is estimated 
as T ^ 10~^s ~ lO^Tosc, while the oscillation frequency 
Tosc ~ lO^^^s is in the THz regime. Although an adi- 
abatic switching of the external magnetic field appears 
necessary, a certain number of quantum manipulations 
seems to be possible. While in metals or semiconductors 
the decoherence will be prohibitively high, defect free in- 
sulators might be much better than this estimate. On the 
other hand, a solid state realization of a memory unit for 
a quantum computer has obvious advantages in terms of 
scalability to devices of higher complexity in comparison 
with EIT based systems. 

(7) While on the one hand the above results are ap- 
plicable to a wide variety of systems, strictly speaking 
the use of the ABC Eq. ( |l5| ) can only be justified in a 
microscopic model, where also the parameter ^ has to be 
determined. This will be accomplished in the following 
for the JPR, because there the problem can be formulated 
as a set of linear finite difference equations and therefore 
a complete solution for all wave vectors can be obtained. 

Thereby it turns out that the optical properties of crys- 
tals with several atoms in the unit cell cannot be de- 
scribed by the function ed^, kz) alone. Then the above 
macroscopic approach based on the slowly varying polar- 
ization PzOi), which is reflected in the ABC Eq. (p^), 
breaks down for both oblique incidence and incidence 
parallel to the layers, see below Sect. IV and V. 



III. MICROSCOPIC APPROACH FOR JPR IN 
CRYSTALS WITH IDENTICAL JUNCTIONS 

A. General equations 

Considering a stack of identical Josephson junctions, 
we label the layers by the index to, the interlayer spacing 
is s and the intrinsic Josephson junctions are character- 
ized by the critical current density Jq. Thus the plasma 
frequency at zero wave vector is given as 



Stt^csJo 
ecO$o 



A? 



(31) 



where $0 is the flux quantum and Ac is the penetration 
length along the c-axis [|lO|-p^. 

In order to determine the transmissivity in the micro- 
scopic approach, we solve the Maxwell equations inside 
the crystal by accounting for supercurrents inside the 2D 
layers at z = ms and interlayer Josephson and quasipar- 
ticle currents, which are driven by the difference Vm,m+i 
of the electrochemical potentials in neighboring layers: 
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N 



771 — 



^dBy 

dz 

Z C 



ms 



(32) 



(m+l)s 



/ (33) 



TV 



ni—0 



(34) 



Pm — Vm,m+1 — esEz^m,m+l + fJ'm+1 — fJ-m- (35) 



Thereby = c/Aab^eao is the in-plane plasma fre- 
quency, Eao is the high frequency in-plane dielectric 
constant and the function / is defined as fm{z) = 1 
at ms < z < (m + l)s and zero outside this inter- 
val. It is seen from Eq. ( p^ ) that the discrete quan- 
tity Pjn — (1/s) jj^^^^^''" Pz{z)dz plays the role of the 
2;-axis polarization Pz{z) averaged between the layers ni 
and m -|- 1, as it describes the response of the Josephson 
plasma oscillations to the electric field in junction m. For 
small amplitude oscillations the supercurrent density is 
given by the phase difference ipm,m+i = 2«eV^_m+i/^ 



as J, 



is) 

m,m+l 



Jo sin (p. 



m,m+l 



m,m-t-l ; 



which was 



used to derive Eq. (34). The difference /im — /^m+i = 
(AiTsa/ ecQ){pm — Pm+i), of the chemical potentials /x,„ 
can be expressed by the 2D charge densities, pm, which in 
turn are related to the electric fields Ez{z — ms±0) near 
the layers by the Poisson equation, ^irpm = Ez{ms + Q) — 



Ez{ms — 0). Further, 



w^(l-i47rCTcw/w^eco)" 



con- 
tains the dissipation due to quasiparticle tunneling cur- 
rents, J^^^+i — fcKn,m+i/es, which are determined by 
the conductivity Cc and driven by the difference Vm.m+i 
of the electrochemical potentials. Note that the assump- 
tion in [ pO| that the quasiparticle current is driven by the 
averaged electric field i?z,m.m+i is an inconsistent treat- 
ment of the dissipation |3^ ]. 

For 2D free electrons we get dp,/ dp = irti^ /{enie) and 
we can estimate the order of a = (eco/47res)(c)/i/9p) 
as ~ 0.38, assuming s = 6.3 A and eco ~ 20. This 
agrees well with a ~ 0.4, which was extracted in the 
one-layer compound SmLai_a;Sra;Cu04_5 from the mag- 
netic field dependence of the plasma peaks in the loss 
function in parallel incidence both in the liquid p8|| and 
the solid phase ||5^. The apparent free electron value of 
the electronic compressibility of the Cu02-layers is not 
in a contradiction to the slightly enhanced effective mass 
m* seen in ARPES measurements |5^, as both quan- 
tities are renormalized differently by interactions. For 
systems with Cu02 multilayers smaller values for the 
compressibility are anticipated due the enhanced den- 
sity of states, effective mass m*, lattice constant s and 
the smaller background dielectric constant EcO, namely 
a ~ 0.05 - 0.1 for Bi-2212 or Tl-2212 (assuming e^o ~ 10 
and d « 12 A), but this quantity can only be extracted 



reliably from experiment. The modification of the dis- 
persion due to nonequilibrium effects is not considered 
in the following, e.g. it is assumed that all frequencies 
are smaller than the charge imbalance and energy relax- 
ation rates [p9|,pU5^. 



B. Dispersion relation 

We obtain now the dispersion relation for eigenmodes 
inside the bulk crystal. For an infinite number of junc- 
tions we average Eqs. ( |3^ ) - ( |35| ) between the layers 
m and m -\- 1 and neglect the discrete layered struc- 
ture, when treating the derivatives with respect to z 
in the Eqs. ( |3^ ) and (|3|), i.e. we replace Ex{z — ms) 

by Ex,m,ni+i = jj;^'^^^^ dzE^ and By{z = ms) by 



B 



y,m-i,m = c?2-By. Usiug the Fourier represen- 

tation with respect to the discrete variable m this gives 
Eq. d) with 



€c{uJ, q) = eco[l - ujl{q)/u}'^], 
ujl{q)^uj^^o[l + 2a{l-cosq)], 



(36) 



(37) 



where < q < 2tt and uj'^{q) describes the dispersion 
of the plasma mode propagating along the c-axis. Using 
Eq. (^) with fc^ = 2(1 — cos q)/s'^, which reflects the ex- 
istence of an upper edge of the plasma band, we obtain 
the dispersion of eigenmodes propagating inside the crys- 
tal in an arbitrary direction. Due to uj'^ea{uj) « ~c^/^ab 
at w « LUcO ^ '^ao we get in the absence of dissipation 
(fTe = 0): 



uj'^{K,q) 



1 



Xlkl 



l + (2Ays2)(l-cosg) 
-|-2a(l — cosg). 



(38) 



The first term in the right hand side of Eq. (^8|) is due to 
the inductive coupling of the in-plane currents excited by 
the component E^ of the electric field. The second term 
reflects the c-axis dispersion due to the charge coupling 
of the intrinsic junctions, which is mediated by variations 
of the electrochemical potential on the layers. For a = 
this dispersion of the plasma mode has already been 
calculated in Ref. |ll| ]. 

For the geometry shown in Fig. |l](a) we can express 
kx = to smO /c via the frequency oj and the angle 9 of the 
incident wave and obtain the dispersion relation for the 
eigenmodes, which are excited by external electromag- 
netic waves. 



= 1 + 2a(l — cosg) 



(a-l)/3 
(5 + 1 — cos q ' 



(39) 



Here /3 = / (2A^^a) ~ 10 describes the inductive cou- 
pling and = 1 — c^A:^/(w^eco) = 1 — sin^6'/eco- To 
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include dissipation, one has to replace uj and w = uj^ juj^Q 
by Cj and w in Eqs. (|3|) and Q. 

In Fig. ^ we plot schematically the dispersion w ~ 
w^/w^Q versus = sin^(q/2). Thereby, z/ is a normalized 
form of the refraction index n = (2c/lus)i' and can be 
used to present both propagating {q real, z/^ £ [0, 1]) and 
decaying {lin{q) ^ 0, i'^ ^ [0, 1]) modes. 



Photon 




-P/2 



^ sin^(q/2) 



Photon 




^ sin^(q/2) 



FIG. 6. Schematic picture of the dispersion relation, 
w — Lu'^/ujco, depending on i/^ = sin^(g/2) (solid line), 
Eq. @, for a = (above) and q 7^ (below). < < 1 
corresponds to propagating solutions with real q, while out- 
side this interval q is complex and the modes decay. It is seen 
that the mixing of a decaying electromagnetic wave (dashed 
line at ly^ = — /3/2) with the plasma band with normal dis- 
persion Q 7^ (dashed) leads to an extremal point Wc and 
a region We < w < a{9), where two propagating eigenmodes 
with normal and anomalous dispersion exist. 

In the absence of charge coupling, a = 0, the eigen- 
mode, which is excited in oblique incidence {a[6) ^1), 
has anomalous dispersion, dLo[v)/dv < 0, cf. Fig. g 
above. It is seen that at a = the width a — 1 of the 
transmission window w € [l,a], where modes can prop- 
agate into the crystal, are determined by the extremal 
values at sm{q/2) = and sm{q/2) ~ 1. 

For normal incidence 9 = a = I) the longitudinal 
plasma mode with a 7^ is decoupled from the trans- 
verse electromagnetic wave as shown by the dashed lines 
in Fig. ^ below, because the electromagnetic wave does 
not have an Ez component which excites plasma oscilla- 
tions between the layers. In this case the wave vector of 
the pure electromagnetic wave inside the crystal is given 
by the relation 1 — cos q + f3 = 0, i.e. the electromagnetic 



wave decays on the scale Xab due to the screening in the 
conducting layers. On the other hand, the wave vector 
of the propagating longitudinal plasma mode, q, is given 
by the relation w = 1 + 2a(l — cos q) and it is real in the 
frequency interval pure 
plasma mode has a normal dispersion, d(jj{v)/di' > 0. 

As a > 1 (<^ ^ 0) is close to unity for any angle 9 
and Eco « 10, the parameter /3 = s^/2aA^j ^ 10~^ ^ 1 is 
small and the two modes mix only when the second and 
third term in Eq. ( p9| ) are approximately equal. This 
happens at small 1/ = sin^(g/2) « u/8a, where the 
small scale u is given as m = [8(a— 1)/3q;]^/^ {u ^ 10^^ for 
cuprates). For any angle 6* 7^ the modes inside the crys- 
tal are a mixture of the longitudinal plasma oscillation 
and the transverse electromagnetic waves. As a conse- 
quence, the electric and magnetic fields of the eigenmode 
are not polarized parallel or perpendicular to the wave 
vector, i.e. the eigenmodes are neither purely transverse 
nor longitudinal. 

From Fig. ^ below it is clear that the mixing of these 
two degrees of freedom at a 7^ 1 and nonzero a can lead 
to the existence of an extremal point We , where the char- 
acter of the dispersion changes and the group velocity 
vanishes. This happens at = 1 + u, provided that 
a > (a— 1)/3/8 and the dissipation is weak, i.e. Im(np) <^ 
Re(rip) or equivalently a = Airac/ujco^co u. We esti- 
mate (7 ~ u in Bi-2212 [ pTf , ct > u in SmLai_a;Srj;Cu04_5 
p8| , p4[ or other cuprates with d-wave order parame- 
ter. Layered s-wave superconductors with the JPR fre- 
quency in the optical interval would be perfect candi- 
dates to study the effects of spatial dispersion, because 
their quasiparticle conductivity is very low at low tem- 
peratures (such systems are possibly realized in organic 
superconductors or intercalated LaSe(NbSe2) [ p^ , 
which has a large anisotropy Bc2.ab/ Bc2,c ~ 50— 130 and 
is therefore expected to be a Josephson coupled system 

In coincidence with the general picture presented in 
Fig. ^ the extremal point uj^ appears near the plasma 
frequency {w = 1), where the wavevector in c-direction 
in the dispersionless theory gets large, see Fig. ^(a). This 
point corresponds to a zero in the dielectric function 
eci^)i as expected from the one mode Fresnel theory, 
cf. Eq. (I). 
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FIG. 7. Real part of ui (above) and 1/2 (below) as a func- 
tion of the normalized squared frequency w/we near the 
plasma resonance for different a {a — Ana /eouJcO,! = 0.26, 
/3 = 10~*, a{9) — 1.1). For We < w < a causality requires 
that Re(i^i) > (Re(i/2) < 0) for the solutions with normal 
(anomalous) dispersion. In the interval 1 — u < w < Wc we 
have in particular Re(i^i) ~ —Ke{i'2), e^ standing waves due 
to interference of vi and U2, see Eqs. MM and (M). 
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FIG. 8. Imaginary part of 1/1,2 for different values of a 
{a — 0.26, /3 = 10"'', a = 1.1). In the region 1 — u < w < We 
below the plasma band we have lm{vi) « Im(;/2). 



In the general case of nonzero dissipation Eq. ( p9D has 
four complex solutions for 1^1.2 at given w = cj^/cj^q, 

vl^iuj) ^ {w - I - 2a f3) /8a ± (40) 
[(w - 1 - 2a/5)2 + 8aP{w - a)]^^'^/8a. 

Near the lower band edge (oj w Wg) this simplifies to 



lyf^iio) = [w - 1 ± y/iw-l)^-u^]/8c 



(41) 



Therefore we obtain 



= Inmal = X%Qu/2as'^ - Xl/{s\ab) > 1 (42) 

in the case of JPR. 

As discussed in Sect. [I A, in a semi-infinite crystal only 
those modes are physical, which decay inside the crystal, 
i.e. Im(i^) > 0, see Fig. ||. For propagating modes this 
implies that the group velocity obeys causality, Vg^ > 0, 
and Re(i'i) > (Re(i/2) < 0) for branches with normal 
(anomalous) dispersion, see Fig. |^. 

We discuss first the limiting case with vanishing dissi- 
pation ((Tc — *■ 0), where the solutions inside the crystal 
are either exponentially decaying (g imaginary) or prop- 
agating modes {q real). For a = we obtain propagating 
mode with real q in the frequency range 1 < w < a 
(cf. dispersion in Fig. |^), and exactly in this inter- 
val the reflection coefficient R < 1. For finite a 7^ 
two physical solutions with real q exist in the interval 
Wc^l + u<w<a provided that a > {a — l)/3/8. In 
the range a < w < 1 -I- 4q! one wave vector, qi, is real 
while the other, i\q2\, is imaginary. The important point 
is that this evanescent solution has small 1921 < 2/3 <C 1 
and because of this it affects strongly the optical proper- 
ties, which are sensitive to large length scales. Outside 
of the interval [we, 1 + 4q!] both gp(w) are imaginary. 

While in the absence of dissipation within the plasma 
band We < w < 1 + Aa at least one of the eigenmodes 
propagates into the crystal, for w ^ 1 we obtain i/^ 2 < 
and the modes qi and 52 decay rapidly on the scales ^/as 
and Xab respectively. 

In the intermediate regime, 1 — m < w < Wg, we have 



Re(gi) = -Rc(g2) = [{u + w- l)/4a]^^^, 
Im(9i) ^ Im(92) = [iu-w + l)/ia]^/^, 



(43) 
(44) 



and the real and imaginary parts of the wave vector q 
are of the same order ^/u (cf. Figs. ^ and ^). There- 
fore, they penetrate deep into the crystal and form 
standing waves, which decay and oscillate on the scale 
[2XabSy/€coo-/ sin 6']^/^. In fact, they are intermediate be- 
tween modes at w <C 1 — u, which decay much faster, and 
propagating modes at 1 + > w > We- 
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C. Eigenmodes of a semi-infinite crystal 

The averaged Maxwell equations (§)-(^) are sufficient 
to determine the bulk dispersion relation Eq. @ of the 
excited eigenmodes and to identify possible critical fre- 
quencies We or uji, where the amplitudes of the excited 
modes equal. At these points the group velocity is ex- 
pected to be low and the microscopic layered structure 
has to be considered more accurately in order to describe 
optical properties. 

For this purpose we solve the electrodynamic equations 
between the layers m and m-l- 1 by using Eqs. ( p^ ) - (35), 
namely the equation 
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By ^ asm6P,n, .g = — ( — ) 
c V a / 



1/2 



(45) 



Physically Eq. (^5|) describes the excitation of a propa- 
gating intrajunction mode with the polarization of the 
electric field in x-direction. Thus at ms < z < (m + l)s 
the solutions for the fields are 

By{z) — Cm e-Kp{igz) + D,n exp(— igx) + asin6'P,„, (46) 
Ex{z) = {eaaay^^^[C„i e^p{igz) - exp(-«gz)], 
-B^(z) = {sin 9 / eco)[Cm e^piig z) + Dm e^p{-igz)] + aPjn- 

Directly from the Maxwell equations follow the continu- 
ity relations, 



Ex{z — ms + 0) = Ex{z = ms — 0), 



By{z — ms + 0) = By{z 



0) + insJx. 



(47) 
(48) 



for the fields By and Ex at layer m with a parallel 
current 'iirJx.m — i^^Ex^z — ms)/uj. Together with 
Eqs. ( ^ ) this leads to the following set of equations 
for Cm = CmG'xp{'i9d{m -\- 1/2)) sin(?/eco and dm — 
Dm exp{—igd{m + 1/2)) sinO /ecQ inside the crystal (A^ — 
2 > TO > 1, A^ is the number of junctions): 

CmV^^ - dmV - Cm-lV + dm-lV^^ = 0, (49) 
2(Cm?7^"^ - Crn-lV) + {a - l){Pm - Prn-l) + 
l{f3/b){CmV'^ -dmV)=0, (50) 
Pm{w -a)+ a{Pm+l + Pm-1 - 2Pm) = 

(sin(fe)/6)(l-2a/3)( ), (51) 

where rj = exp(i6) and the small parameter b = gs/2 
s/Xc ^ 10~^ ^ 1 characterizes the discreteness of the 
crystal structure. We will assume in the following that 
/3, & < and g - P^/^ - 2/Aafc, as it is fulfilled for 
highly anisotropic (Ac ^ Xab) layered superconductors, 
e.g. Bi- or Tl-based cuprates. In our calculations we will 
keep only the terms of lowest order in the small parame- 
ters l3 and b. Eqs. (H9) - (M) give the dispersion relation 



Eq. (39) with high accuracy b^ and a/3. This difference 
between the exact result following from the Eqs. ( ^ ) to 
(O) and the averaged dispersion (cf. Eqs. (p6|)) can be 



understood explicitly from Eqs. (pq): the replacement of 
Ex{ms) by the averaged Ex^m,m+i is correct in order 6., 
i.e. when neglecting the discrete layered structure within 
the unit cell. 

The solution inside the crystal has the form 



^ 7pexp(igpTO), 

p=l,2 



E 

= 1.2 



dm. = ' jpd{qp)e^p{iqpm), 



Pr 

d{q) = 



lpP{qp)e^p{iqpm), 



p=l,2 

\ — rf' exp(— ig) 
if- — exp(— «g) ' 

1 + d{q) 
w ^ a — 2a{l — cos g) ' 



(52) 
(53) 
(54) 

(55) 
(56) 



where qp{u;) are the wave vectors of the eigenmodes for a 
given frequency as determined by Eq. (^ ) and 7p denote 
the relative amplitude of the excited modes, which is to 
be determined next. 

Neglecting the layered structure, e.g. b ^ s/Xc 
and — > 1, we obtain c = d = 1. In this case we 
can relate the variables c™, dm and Pm with the elec- 
tric and magnetic field averaged between the layers, i.e. 
Ex,m,m+i oc Cm + dm and Ez^m,m+i is mainly determined 
by the polarisation Pm- 



D. Microscopic Boundary condition 

Now we find the ratio of the amplitudes 71 and 72 
microscopically by solving the electrodynamics of the 
surface junctions explicitly rather than using any phe- 
nomenological ABC. The equations for the first super- 
conducting layer (to = 0), which arc complementary to 
the Eqs. (||)-(|53), read as 



co + do + {a- l)Po = 
cor/^^ — dor/ sin 29 



sm9 



B 



rcf\ 



(By - Bf), 



Po{w -a) + a(Pi - Po - aP„) - 

(1 + a)(co + do) = asin0(P;" + Bf). 



(57) 
(58) 

(59) 



y , Bf are the magnetic fields for incident and 



Here B, 

refiected light respectively. We omitted in these equa- 
tions terms proportional to /3/sin(g/2) and 6/ sin(g/2), 
which are of order /3^/'' ^ ^a/n-p ^ (s/Xab)^^'^ *C 1 and 
b/13^/'^ ~ 1/n? ~ sXab/X^ ^ 1 in comparison with re- 
maining terms of order unity. After eliminating the fields 



Byi and By^ from Eqs. (|5^) - ( [59| ) we obtain in lowest or- 
der in b and [3 the microscopic boundary condition 
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Po{w - a) + a{Pi - 2Po) - (60) 
a(eeo + l)[co + do + {a- l)Po] - {cq + do) = 0. 

We can present this condition in a more transparent 
form by calculating the difference between Eq. ( Jsi] ) for 
771 = 0, where P_i (outside the crystal) is formally given 
by Eq. ( [5^ ) for 777 — —1, and the real equation for Pq, 
Eq. (|60|). We also take into account that in the lowest 
order in /3 and b we obtain the relation cg+do + aPo ~ Pq 
near We with accuracy ^ (ea/77p)^^'^ ^ 1 using 

Eqs. ( |55| ) and ([s^). This gives the boundary condition 

Pm=~i = ^ ^pP{qp)e-Kp{-iqp) = 0, (61) 

p=l,2 

which has the simple interpretation that the surface junc- 
tion (m = 0) has only one neighboring junction, i.e. the 
junction ?n = — 1 is absent. This result is a microscopic 
derivation of the ABC Eq. (15) by noting that Pm is the 
average macroscopic polarization Pz{z) between neigh- 
boring layers, i.e. 



P-^ = 



1 



P,(z) «P,(z = 0)-s9,P,(z = 0). (62) 



Taking into account that the deviation of R from unity 
is significant only when \qp\ <C 1, we expand Eq. (61) in 
Qp by using P{q) — P(0) ~ q"^ (in leading order in b) from 
Eq. (p6|) and obtain Eq. (|l6|) with £ = —s: 



p=l,2 



«9p) = X! 77)(l + *C"p) = 0. (63) 

p=l,2 



Note that this result and consequently also the expression 
for K, Eq. (17), is only valid in leading order in ea/rip ^ 
« 1. 

With this identification of the parameter ^ we can es- 
timate 



^|77,i772| Ac/Aah > C^a ^ sK/>ib > 1 



(64) 



at to = LUe in Bi- and Tl-based layered superconduc- 
tors. This shows that when the anisotropy Ac/ A^fc is large 
enough, the atomic structure modifies strongly the trans- 
mission, cf. Eqs. (^ and (|l^). Here we also justify the 
relations discussed in Sect. 11, 



i\nin2\ < 771 



A. 



n2 



(65) 



which allow us to neglect the atomic structure away from 
a small frequency intervall of width ^ u^^^ around iVe. 

Due to |77i| > |772| (cf. Fig. 0) and |7ii772|e^i ~ 
^ab/ 3 ^ 1 away from uOf, or uoi the usual one- mode Fresnel 
theory is valid everywhere, except near the resonances at 



E. The transmission coefficient 

As a consequence, we reproduce in our microscopic the- 
ory Eq. ( |l7| ) for K and therefore the transmission and re- 
flection coefficients T, Tp, ppp', Eqs. ([l9| ) and (|2^ ) - ([2^). 
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FIG. 9. The dependence of the real and imaginary part 
of Ko ~ k(cos 6'eaC/2) on w = uj'^/ll>1(, for cr = 5 10~* and 
a = 0.001. The lineshape of tcq is asymmetric with a sharp 
edge at the extremal frequency lOc = 1 -I- w « 1 and the upper 
edge at w = a{6) = 1.05, which is determined by the angle 6 
of incidence. 
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FIG. 10. Transmission T depending on w/we near the 
JPR frequency for conductivities a = lO"'^ (left) or cr = 0.01 
(right) for various a {f3 — 10~*, a — 1.1, ^eaC0s9 = 2). For 
low dissipation a <^ (left) the resonance is additionally 
damped due to a in the region near We, where the Fresnel 
approach is invalid. 

The real and imaginary parts of k = e^g are shown in 
Fig. 1^ and have a characteristic shape with a sharp edge 
at the extremal point We, provided that the dissipation 
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a is small. The real part is dominant only in the interval 
We < w < a, where both modes vi and 1^2 are real and 
propagating and the transmission T into the crystal is 
significant. The window of transmission We max ~ — 
w < a is therefore only determined by a{9), and not by 
the bandwidth ^ a. The width of the peak in k and 
T near Wc,inax assuming a 10% criterion is of the order 
lOOu. 

In the interval a < w < 1 + 4a, where vi ~ 
becomes imaginary and small, while 1^2 is real, k is a 
complex number (even for ac = 0) with a real part pro- 
portional to ji^^'l. In contrast to the standard Fresnel 
expressions, this makes transmission possible, but it is 
weak of the order 6, because only a small part of the in- 
cident light transforms into a propagating mode. There- 
fore, deviations of R from unity are significant only in 
the frequency range We ~ 1 ^ w < a, as in the system 
without dispersion. 

If the dissipation is very weak. 



Im(ni + n2) <^ ^"1^2 4^ a ^ u'^ 



(66) 



the nonuniversal term characterized by the parameter ^ 
in Eq. ^ is important. Then according to Eq. ( p^ the 
maximum of T is reached at Wcmax = + u^^'^ / VSa. 
The amplitude, 



IV. CRYSTAL WITH ALTERNATING 
JOSEPHSON JUNCTIONS 

For the geometry in Fig. |l|(a) we consider the crys- 
tal with two alternating Josephson junctions I — 1,2 
characterized by different critical current densities Jq.i 
and two bare plasma frequencies Wco.i and Wco,2 related 
to Jo,i and Jo,2 as described by Eq. (|3|). We de- 
note w — uj'^/uj'^Q i and S ~ ^c0,i/^c0,2 < 1- the 
view of recent experiments ]3^ , we also allow for dif- 
ferent c-axis conductivities ct; {I — 1,2), which are ex- 
pected to vary according to the different tunnel ma- 
trix elements in the junctions, ai/a2 — '^c0i/'^c0 2i as 
found for La2-a;Sra;Cu04 |Q, and which are assumed 
to be frequency independent in the following (ct; = 
47rcr//ecoWco,i). AH other parameters of the junctions are 
assumed to be identical. 

The equations inside the crystal are analogous to 
Eqs. (^)- (|5l|) and the details of their solution are given 
in appendix |A[ Here we summarize the main features on 
the basis of the schematic dispersion w{v'^ = sm^{q/2)) 
in Fig. |l^ and the squared refraction indices lyf 2{w) in 
Fig. @. 



T 



[l + (sAc COS 0/(V^A2J)2] 1/2 + 1' 



(67) 



is smaller than unity and it depends on the microscopic 
structure via the factor sAc^/cdj/A^^ which may be of or- 
der unity in cuprates like Tl-2212 with Xc/Xab ^ 100 and 
the JPR frequency 20 cm^^. This effect can be seen 
in Fig. |lfl (left): Without dispersion, i.e. a '3> u{a)'^ for 
a — 1(P^, the peak amplitude is limited by the small 
dissipation, a, only, while for a = 0.1 {a ^ u^) the 
peak at Ci^cmax is damped additionally due to the novel 
term i^nin2 in Eq. (p7|), as discussed above. Physically 
this can be understood from the fact that the vanishing 
group velocity leads to a slow motion of the wave-packet 
and hence makes the transmission sensitive to the in- 
homogeneous layered structure of the system, i.e. the 
translational invariance of the system is broken. 

On the other hand, high dissipation ct 3> m overshad- 
ows the effect of spatial dispersion completely (Fig. [ic| , 
right). In this case the result near the lower edge of the 
transmission window is almost the same as in the disper- 
sionless model. 
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c cos 61 



LuXn 



1/2 



(68) 



and is mainly determined by ac- 




-2P 



^ sin2(q/2) 



FIG. 11. Schematic picture of the dispersion 

w{u^) = uj^{i')/lj^q i for two alternating junctions (ctc — 0). 
The dispersion of the plasma mode at 6' = (dashed line), 
i.e. when it is decoupled from the electromagnetic wave, is 
normal in the lower band. Its mixing with a decaying electro- 
magnetic wave (as shown by the dashed, vertical line at nega- 
tive = sin^(g/2) ^ —213) results in two propagating modes 
(solid) near the lower band edge We- This frequency forms 
an extremal point with vanishing group velocity as in the one 
band case (cf. Fig. ^ below). The anomalous dispersion in 
the uncoupled upper band gives rise to one propagating and 
one decaying mode and a special point Wi, where qi = —iq2- 



The band edges Wj 



± 

low, up 



are defined in the appendix A. 
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ima we obtain approximately 



We 




(l + (5)(l + 2a)/25T (69) 
[(1 + 5)2(1 + 2a)2 - 45(1 + 4a)]i/V25. 



As for identical layers the optical properties are dom- 
inated by the mode with smaller \np\ = c\vp\/sLU and 
significant deviations from the one mode Fresnel regime 
occurs at w 1 712!. 




We Wq wl^l 



FIG. 12. Schematic dependence of the squared refraction 
indices i/f (solid) and (dashed) on the squared frequency 
w = tj^/tJcO,! for two alternating junctions {ai — > 0) with the 
peak positions We,i and the band edges as defined in 

appendix 

At perpendicular incidence a = 1 [9 = 0) the lon- 
gitudinal plasma mode is decoupled from the transverse 
electromagnetic wave, as the incident electric field has 
no component perpendicular to the layers. In this case 
the lower (upper) plasma bands have a normal (anoma- 
lous) c-axis dispersion (dashed lines in Fig. |l^) due to 
the charge coupling a. 

In contrast to this, for a > 1 the frequency uj{q) in- 
creases as g — > due to the inductive coupling in both 
bands (solid lines in Fig. |ll|). For the lower band this 
can lead for sufficiently large a to an extremal point We 
at the lower band edge as in the case of identical layers, 
where (for cr/ = 0) two modes with real q exist, while near 
the upper band edge w^^^ o^i^ mode propagates and the 
other decays. In the upper band there is one real and 
one imaginary solution everywhere in the band due to 
the anomalous dispersion, and according to Eq. ( pO| ) in 
the general section II the maximal transmission is at Wi , 
where ~ —q^- 

All special frequencies mentioned in Figs. |ll| and |l^ 
are explicitly expressed by microscopic parameters in ap- 
pendix ^ For the frequencies w^^i of the resonance max- 
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FIG. 13. The real part of the effective dielectric func- 
tion Ecft ~ without explicit dispersion (a = 0, 
(5 = t^co,i/^co,2 = 0.3, a = 1.1, /3 = 10"*, J; = 0). In tliis 
limit the dielectric function is directly related to the refrac- 
tion index and the wavevector of the single excited mode, 
Eeff C)C no oc fc^, see Eqs. (^ and (^. The poles at the lower 
band edges, where the averaged e^iy^) — vanishes, indicate 
the appearance of the special points We and Wi in the two 
mode theory, cf. Fig. |l2| for the corresponding case a 7^ and 
Fig. ^ for the general picture. 

Keeping only the solutions with smallest |i'p| near Wi 
and We, e.g. 1^2 (1^1) for w < Wi {w > Wi) in the upper 
band, we obtain in the limit a — > a pole in ^ q^ ^ 
Ecff, as can be seen from Fig. |l^. This is an explicit 
microscopic confirmation of the general expectation that 
critical frequencies, where \v\ \ = \v2\, appear, if there is a 
pole in fc^ ^ Eoff. For oblique incidence the singularities 
in ecff(w) coincide with the zeros of the averaged e(w) 
introduced in Eq. (||) with td = eco(l — Joj"^) ■ This 
can be expected from a macroscopic treatment, where 
the spatially averaged e is introduced in Eq. 

The regularization of the poles in eoff is seen by com- 
paring the behaviour of Re(fp) ~ Re{klp) for a ^ 
(Fig. |l2|) and of Re(eoff) - fc^ for a = (Fig. |l3|) near 
We and Wi with the schematic picture in Fig. |^. 

In sect. V it will be shown that a situation, where 
a second mode contributes in a similar way as near the 
point Wi in the upper band, can also develop from a pole 
in the dispersionless dielectric function without explicit 
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spatial dispersion, e.g. for a = 0, due to the intrinsic 



atomic structure within the unit ceU, see Fig. 16. 

Also like in the single junction case, the transmission 
into the crystal in the lower band is only significant, 
if both excited modes are propagating into the crystal. 
Consequently, the width of the resonance Aiow ~ "i^o ~ ''^e 



(Aup 



i) in transmission T(w) in the lower (up- 



per) band are considerably smaller than the band width 
of the allowed eigenmodes in the crystal. 



low 



Wp or 



w+p — w~p respectively, see Fig. |l^. 

As derived in Appendix ^ the additional boundary 
condition near the special points We and Wi is analogous 
to the case of identical junctions Eq. (^l]) and reflects the 
fact that on the surface one neighboring junction is miss- 
ing. In leading order of (3^^^ ^ ^a/n^ and 5//?^/'' ~ 1/n^ 
we obtain 

-Pm=-i,2 = ^ 7pP2(gp) exp(-igp) = 0. (70) 

p=l,2 

Thereby Pm=-i,2 = /^^ Pz2dz is the average of the I — 2 
component Pz2 of the macroscopic polarization vector in 
the missing junction in the cell m — —1, {Pi{q), P2{q)) 
denotes the eigenvector of the excited mode and jp 
describes the rel ative amplitude of the excited modes 
p = 1,2, see Eq. ( A13 ). This microscopic result gives an 
a posteriori justification of the phenomenological ABC 
in Eq. ( |l5| ) for the multimodc case, where the length 
scale £ — —2s is identified with the lattice constant in c- 
direction. This shows in particular that the macroscopic 
approach is possible, if and only if multicomponent local 
polarizations Piz,P2z inside the unit cell arc introduced. 
Expanding ^(g) ^2(9) exp(-ig) » P2(0)(l - iq/2) 
with the help of Eq. ( [A1^ ) for ^2(9) ~ ^2(0)(1 -I- iq) and 
taking into account the doubled unit cell in g = 2skz , we 
obtain Eq. ( p^ with the effective parameter ^ = —los/c. 
It is pointed out that the same result for the amplitude 
ratio of the excited modes as in the single layer case is 
reached here in a nontrivial way by an interplay of the 
lattice constant £ = 2s and the internal structure of the 
eigenmodes contained in Pi{q). 




FIG. 14. Real and imaginary part of = k(cos 0eaCJs/c) 
for alter- 
nating junctions (a = 0.2, a — 1.1, 5 = (j-'c0,i/'^c0,2 = 0.3, 
P = 10""*) near ~ 1.244 (left) and Wi ^ 4.823 (right) for 
different quasiparticle dissipation ai = So-2 = 10"^ (above) 
or 0.005 (below). 

Now we are in the position to calculate the reflection 
and transmission coefficients R and T near the resonances 
in leading order in f]^/^ ~ ^a/n^ and hjfi^l^ ^ I/'t^^, 
where k is given by 



nin2[P2('7i)»^2 - P2{<l2)ni\ 



eaCOS0 P2{qi)nl~ P2{q2)nl 



(71) 



and shown in Fig. |Tj. Due to the lattice constant 2s the 
refraction indices of the bulk eigenmodes are here rip — 
ckz.p/to — cVp/sLJ. This result reduces to Eq. ^ when 
expanding P2{q), and the results of the general Sect. || 
can be used. 




1 1.01 1.02 1.03 0.995 1.005 1.015 1.025 



0.08 
T 




0.95 



0.06 



1 1.05 




17 



FIG. 15. Transmissivity T near the lower (left) and upper 
(right) plasma bands with the frequency axis normalized to 
wiow (wup) respectively. Parameters: S — 0.3, [3 = 10"^, 
ea^cos6' = 1, different conductivities in the plots above 

5^2 = 0.01, 



(ai = Sa2 — 10 a = 1.05) and below (cti 
a = 1.1) and varying a (see plot). 



The lower band is similar to the case of identical 
junctions in the sense that in the transmission win- 
dow We < w < wq two propagating modes are excited 
and we obtain the same maximal transmission coefficient 



max, low 



(cf. Eq. dl 



In the upper band we get from the general Eq. (20) for 



small dissipation ai ^ u^p at w — Wi 



- max, up 



cost 



(g - l)£up 
8a^a 



1/4 



(72) 



which is smaller by the factor {s/\ab)^^^ than Tmax.iow 
(-^up = 'Wi{l + S) — 2 — 8a). This can be seen in the 
upper part of Fig. |l5|, where for low ai <C wiow,up ^ 1 
(see definitions in App. ^) the upper plasma resonance is 
considerably suppressed by increasing a, while the lower 
band is weakly affected. This suppression can be un- 
derstood physically by the fact that at the surface the 
energy of the incident wave is distributed between a prop- 
agating wave and a decaying (and finally reflected) one 
and is therefore less efficiently transmitted in the crystal 
than in the lower band, where the two excited modes are 
propagating. Physically, the eigenvectors near q « in 
the lower (upper) band involve in phase (out of phase) 
plasma oscillations and consequently external long wave 
length radiation couples more efficiently to the excita- 
tions in the lower than those in the upper band. 

The difference between the values of Tmaxjiow and 
Tmax^up decreases as dissipation increases, see Fig. |l^ 
below. It vanishes in the Fresnel limit, for which 
{Airai / uJcO,i^co){^ab/ C(s) becomes much larger than unity. 

In oblique incidence the suppression of the peak in the 
upper band is quite limited to systems with very low dis- 
sipation and perfect crystal structure and might be diffi- 
cult to observe in SmLai_2,Sra;Cu04_5. Instead of this, 
a quite high ratio of the peak amplitudes has been ob- 
served in this material for incidence parallel to the layers 
pll-pSl, see below and Ref. W^M. 



V. INCIDENCE OF LIGHT PARALLEL TO THE 
LAYERS 

In this Section we discuss the reflectivity for incidence 
parallel to the layers, cf. Fig. ^b) for 9 = 0, in the 
crystal with two alternating junctions, when the explicit 
spatial dispersion, i.e. the dependence on the wavevector 
q, is negligible. 



We will confirm microscopically the breakdown of the 
macroscopic Fresnel approach using the effective dielec- 
tric function Cc, Eq. (^, when the wavevector {k^pl of 
the excited modes becomes large and the group velocity 
is small, cf. Eq. (|l|). This happens near the pole Wpoio 
of ic, which coincides with the upper edge of the lower 
plasma resonance in the reflectivity (cf. Fig |l^). This 
frequency is sometimes associated with the excitation of 
a so called "transverse" mode p6| , p7t , although all the 
modes excited in the plasma bands are transverse in this 
geometry. For simplicity we will present here the formu- 
las for a = 0, the general results are given in Appendix 

i . 

Physically, the conventional theory is insufficient, be- 
cause it averages the Eqs. ( ^2| ) - (35) within the unit 
cell and neglects the electric field components parallel 
to the layers, in order to arrive at the response function 
ic for the averaged field Ez.m,av = Jms^^^'^ E^dz. This 
corresponds to neglecting the average jj^^'^^^'^ dzdzBy = 
By[{m + 2)s] — By{ms) and the average of dzE^ respec- 
tively, i.e. to setting fc^ = in the Eqs. (||) - (||). This as- 
sumption is justified away from Wpoic, where the wavevec- 
tor \kx\ is small, as the gradient of the electric field van- 
ishes, if the charge density on the layers is slowly varying, 
dzEx ~ \kx\- On the other hand, at Wpoio the charge den- 
sity varies on atomic scales, the intra-junction mode with 
polarization of the electric field in x-direction is excited 
strongly and the basic assumption of the averaged theory 
is invalid. 

A more careful averaging of the Eqs. (|3^)-(|35|) within 
the junctions rather than the unit cell leads to relation 
between average electric fields inside junctions {a = 0, 
CT/ = 0, for the general case see App. ^) 



ckx 

UJ 



2 + /3o 



l + /3o 1 
1 l + /?o 



Ez,\ 

Ez,2 



= 0. 



(73) 



Here the dielectric tensor e is given as eu = Cd = eco(l — 
c^2^/w2), = £21 = and (3o = aP = 3^/(2X1,,) < 1 
accounts for the coupling of the averaged electric fields 
Ez.1,2 in the junctions of type Z = 1,2 via the electric 
field component E^- The latter is weak, ~ (3o, due to 
the strong anisotropy of the material. For /3o = one 
eigenmode of Eq. ( [73|) corresponds to the solution in the 
averaged theory determined by c^fc^^ = w^ec, and its 
eigenvector obeys Dz = Cci^-zi — ^c2Ez2 as it is assumed 
in macroscopic electrodynamics. Consequently, near the 
lower band edges ujca,i only the plasmon in the junction 
of type I is excited. The other mode has an eigenvalue 
l/kx2 — and corresponds to an out of phase mode with 
the eigenvector Ezi = —Ez2, which is not excited by a 
homogenous incident beam. Accounting for the excita- 
tion of the electric field components parallel to the layers, 
at /3o > 0, both modes mix and the singularity at Wpoic is 
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removed. This is shown schematically in Fig. and is a 
consequence of the presence of two junctions in the unit 
cell (implicit spatial dispersion), even in the absence of 
c-axis coupling (a = or g = 0). 

Including the dissipation due to ai, in the case q <^ b 
sinO ^ 1) of parallel incidence the dispersion is given 
by Eq. ( |A5| ), and the general solutions rixp — ckxp/uj are 
presented in Eqs. (Bl) in the appendix. Away from the 
pole ■Wpoie these solutions can be expanded in leading or- 
der in /3o and we obtain the usual wave with the refraction 
index Uxi = ckxi/to corresponding to the average ec(^) 
(for a = 0): 



1 



1 



c k^i _ ec{w) _ 5{w - wiow)(w - Wup) + iS 



5i = (1/2)u;3/2j(cti +^2) 



(74) 



S 



w 



1/2 U 



<7i{5w - 1) + 6a2iw - 1)], 



The zeros of fc^i are at the plasma edges wio 



u;iow,up= (l + <5)(l + 2a)/2(5T (75) 
[(1 + 6f{l + 2af - 4^(1 + 4a(l - Pz^))]"'^ /25. 

For /3o = this corresponds to the single excited mode 
fc^ ~ ic and we see from Fig. ^ (dashed line) that \nxi\ 
becomes large at the pole Wpoic — (1/2 + 2a)(l + 5^^)- 
The discrete layered structure (/3o ^ 0) results in the 
regularization of the pole and its transformation into a 
special frequency uji, where = without dissipa- 

tion, see Fig. |l^ (solid) . This is similar to the behavior in 
the upper plasma band in oblique incidence, see Fig. |l^, 
where a pole in the one mode Fresnel dielectric function 
Eoff is transformed into the special point w;. There it was 
a consequence of explicit spatial dispersion (a 7^ 0), while 
now the second solution 71^2 appears due to the atomic 
structure within the unit cell even at a = 0. 
Away from the pole Wpoic we get (for a = 0) 



1 



1 

0-2 



w6{w - Wpoic) + iSi 



(76) 



and \kx2\ ^ 0(1/ Po) is large in comparison with \kxi\. As 
the solution with the smallest refraction index determines 
the optical properties, the wave with wave vector kx2 can 
therefore be neglected everywhere except at Wpoic, where 
kx2 is small at weak dissipation and the general Eqs. (Bl) 
have to be used. 

Let us now find the solutions for the magnetic and 
electric fields inside the crystal which determine the re- 
flection coefficient Rn = |(1 — k\\)/{1 + 




FIG. 16. Schematic refraction index n^(it;) without dis- 
persion Q = 0, but with (solid line /3o ~ s^/2A^5 > 0) or 
without (dashed, /3o = 0) accounting for the intrinsic inho- 
mogeneity in the unit cell. The latter implicit spatial dis- 
persion corresponds to the excitation of the mode with elec- 
tric field polarization parallel to the layers. The frequencies 
'w{kx = 0) = wiow.up form the plasma edges in the reflectivity 
R\\ . We can also interpret the plot for /3o = as the averaged 



dielectric function tc 



kl ~ nl (cf. Eq. (0)) as a func- 
tion of w. Then the pole in edw) at Wpoic in the one mode 
approach indicates the appearance of a special frequency for 
/3o 7^ 0, where n^i = —n^2, which is similar to the general 
picture in Fig. |^ and the upper band in oblique incidence, cf. 
Fig. |l|. 

The solution for By a,t x < consists of the incident 



and reflected waves Bl] 



and Bl''^ 



,y ^ii^ , which are homogeneous 
in the z-direction, and the wave By with k^ 7^ 0, which 
is excited due to the inhomogeneity of the crystal in z- 
direction and which is localized near the surface: 



By{x) = B'^ exp(ifc,x) + Bf exp(-zA:,x) 
-I- ^ By{kz)eyip{ikzZ - ik^x), 



(77) 



where c^(fc^ + A;^) = oj'^. The solution at a; > is given by 
Eq. (^6|), when introducing k^p explicitly by substituting 
sin 6 — > —ckx/io and taking into account the superposi- 
tion of the two solutions p = 1,2. 

In addition to this, we need an additional boundary 
condition, in order to determine the ratio, in which the 
modes p — \,2 are excited. Kt 9 — Q the in-plane currents 
Jx,m = i^1,Ex{z = ms) /[Airuj] and consequently the 
components inside the layers m vanish at x — s- 0. This is 
equivalent to Pekar's boundary condition Pj; = 0, which 
turns out to be sufficient in this case due to the absence 
of extremal points, where ni -I- 712 ~ (cf. Fig. p6[ ). 

As worked out in appendix^ this leads to the reflection 
coefficient R\\, Eq. (^, where (a — 0) 

a2Z 



— /<rCff — 



" ainxi + a2nx2Z' 



(78) 
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z = - 



fp 



(l-/i)(l + /2)(ai-l)fcx2a2 
(l-/2)(l + /i)(a2-l)fcxiai' 
{wi - l)(2ap + /3o) - {ttp - l){ap 



(flp - l){ap + Po) 



(79) 
(80) 
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FIG. 17. Reflectivity = jl - + in paral- 

lel incidence for a = 0, 5 = '^c0,i/'^c0,2 = 0.3, EcO = 19, 
/3o = s^/2A^i, = 10~* as in the cuprates and difi'erent conduc- 
tivities CTi = 55-2 (see plot). 
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FIG. 18. Reflectivity ity = 
incidence for a = 0, 5 = 0.3, a 



|1 -f Kii I in parallel 
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19 and different 
(Ho — s /2Ao6 (see plot). It is seen that, when the second 
solution is taken into account, at /3o 3> 0.0001, the reflectivity 
drops near the pole Mpoie of the averaged dielectric function 
Ec, where n^j — — n^2- 

In Fig. the reflectivity R\\ [w) is shown for different 
conductivities di = Sa2 and a value /3o = 10^"* appro- 
priate for high temperature superconductors. The reso- 
nances in the lower and upper plasma bands are asym- 
metric and have a sharp lower edge at wiow.up- The upper 



edge of the lower band is given by the pole Wpoio, where 
in the conventional one mode theory £c becomes nega- 
tive and the single excited mode with an imaginary wave 
vector decays (cf. Fig. |l^) . 

In Fig. |l^ the effect of the discrete layered structure 
(/3o 7^ 0) on the reflectivity in parallel incidence is 
shown. At the special point Wpoie, where the second so- 
lution — becomes relevant, the reflectivity i?|| 
drops for large Po and the lineshape is modified. This 
behavior is similar to the resonance at Wi in the upper 
plasma band for oblique incidence due to spatial disper- 
sion, see Sect. IV. This modification of the JPR lineshape 
is beyond the conventional one mode Fresnel approach, 
which is valid away from Wpoioi in particular near the 
plasma resonances. 

Therefore for the interpretation of the main peak am- 
plitudes the simplified effective dielectric function Cc is 
sufficient and has been used in to extract the pa- 
rameter a « 0.4 from the experimental loss function in 
SmLai_a;Sra;Cu04_5. In contrast to the dissipation 
was introduced here microscopically in the quasiparticle 
currents and it is taken into account that the quasipar- 



ticle conductivities alternate, (Ti/(72 



cO 



c0,2, in the 



same way as the critical current densities and plasma 
frequencies Wco,;- Correctly accounting for dissipation is 
crucial for a quantitative interpretation of the experimen- 
tal loss function. As the parameter a can be extracted 
independently from the magnetic field dependence of the 
plasma resonances, see [ p8[ , this is also a way to deter- 
mine the c-axis conductivities cr/. Both ways to extract 
a « 0.4 from far-infrared data are well compatible with 
the ARPES measurements [Esf. 



VI. CONCLUSIONS 

In conclusion, the effect of spatial dispersion and the 
atomic structure on the optical properties of stronly 
anisotropic uniaxial crystals has been studied in general, 
taking as a generic example the Josephson Plasma Res- 
onance in stacks of identical or alternating junctions. 

Thereby, multiple eigenmodes, propagating or decay- 
ing, are excited by incident light, which interfere with 
each other. This intrinsic birefringence can be detected 
in transmission by oscillations with respect to the sam- 
ple thickness or the splitting of the incoming (laser) beam 
(cf. Sect. II.C). 

In contrast to the usual assumption that the effect of 
dispersion or of the atomic structure on optical charac- 
teristics is strongly suppressed ~ s/X <C 1, as the wave- 
length A of light is much larger than the lattice constant s, 
we showed that near resonance frequencies the refiectiv- 
ity may differ significantly from the conventional Fresnel 
formulas, if dissipation and disorder are weak. 

Near extremal frequencies LUe, where the group veloc- 
ity Vg = Xg/uj vanishes, the stopping of the wave packet 
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makes the propagating light sensitive to short length 
scales Xg. As a consequence, for oblique incidence the 
transmissivity into the crystal cannot be expressed by 
the bulk dielectric function alone and the amplitude of 
the resonance near LUg crucially depends on the atomic 
structure of the crystal. This additional damping due 
to the c-axis coupling a for low dissipation is shown in 
Fig. In contrast to this, the width of the resonance in 
transmission is not affected by the c-axis charge coupling, 
but is rather determined by the angle of incidence. 

These extremal points uje may appear, whenever an 
optically active crystal mode with normal (anomalous) 
dispersion is mixed with a propagating (decaying) elec- 
tromagnetic wave. For these results it was crucially to re- 
alize that the resulting two eigenmodes with normal and 
anomalous dispersion have wave vectors and refraction 
indices with opposite sign near We in order to preserve 
causality. 

In addition to this, for a crystal with several opti- 
cal bands we predict different amplitudes of the reso- 
nance transmission into bands, which are characterized 
by different types of dispersion and which are equivalent 
in a dispersionless theory. When inside the crystal one 
mode is propagating and the other one is decaying, the 
maximum of T is at frequencies uji, where the relation 
ni = —in2 for the refraction indices holds. At these fre- 
quencies the peak amplitude of T is strongly suppressed 
in comparison with bands, where the two excited modes 
are propagating (Fig. |l^), provided that the dissipation 
is low. This provides the unique opportunity to extract 
microscopic information about the eigenvectors of the ex- 
cited modes from the line shape in optical experiments. 

For incidence parallel to alternating layers a second 
mode is excited even without explicit spatial dispersion 
(k-dependence of Cc) due to the intrinsic inhomogeneity 
within the unit cell. Near the pole of the effective di- 
electric function at the upper edge of the lower plasma 
band a special point appears, where n^^ — ^'^^ 
an appropriate choice of parameters this can modify the 
lincshape of the resonance. 

This behavior near uje and uJi cannot be obtained in 
the one mode approach without dispersion. The only 
intrinsic indication for the breakdown of the conventional 
Fresnel theory is the appearance of poles in the effective 
dielectric function ecs, see the schemcatic Fig. |[ There 
the excited wave vectors ~ are large, the group 
velocity is small, cf. Eq. (|l|), and concomitantly small 
atomic length scales become important (cf. Fig. |l^ and 
[TgI for oblique and parallel incidence). 

These features were demonstrated explicitly for the 
JPR with identical and different alternating junctions, 
but they are general for any modes, e.g. for optical 
phonons with anomalous dispersion in insulators, which 
form a polariton branch with an extremal point, see 
Fig. 1^. However, the condition of weak dissipation and a 
perfect crystal structure are crucial to observe deviations 



from the Fresnel regime. 

For the JPR this theory was used to extract the pa- 
rameter a w 0.4 from the optical data obtained for 
SmLai_a:Sra;Cu04_5 with two different alternating in- 
trinsic Josephson junctions between the Cu02 single lay- 
This value corresponds to an electronic com- 
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ers 

pressibility, which is unrenormalized by the interaction, 
while for multilayer cuprates a smaller value of a is ex- 
pected. This result is compatible with the ARPES mea- 
surements [ p5[ and gives an important input parame- 
ter for the coupled Josephson dynamics in the stack. 
Thereby the correct treatment of the c-axis conductiv- 
ities in different junctions is essential for a quantitative 
interpretation. 

It is also pointed out that spatial dispersion provides 
a way to stop light in a crystal, which is dual to previ- 
ous proposals based on the frequency dispersion of the 
medium, see Sect. II. D. From the application point of 
view, this suggest future magneto-optical devices (using 
e.g. the JPR) for storing light coherently, as it is required 
in an optical quantum computer. By imprinting a group 
velocity profile with the help of an inhomogeneous ex- 
ternal magnetic field, event horizons with respect to the 
propagation of light can be created in a solid. 

To summarize, possible experiments to demonstrate 
the effect of spatial dispersion on the optical properties of 
solids include the demonstration of: (a) intrinsic birefrin- 
gence and beam splitting, (b) stopping (delaying) light 
pulses, (c) the relative amplitude of bands with a dif- 
ferent number of propagating excited modes and (d) the 
intrinsic damping of peak amplitudes in materials with 
negligible dissipation and disorder. 

From a general point of view, these results shed new 
light on the long standing question of the treatment of 
spatial dispersion for optical properties of solids and pro- 
vide the first microscopic derivation of the ABC as sug- 
gested in H]. It is expected that the phenomenological 
results presented here can have wide implications for the 
interpretation of resonance amplitudes and lineshapes in 
optical experiments, especially near frequencies Ue or LUi, 
which appear near poles of the conventional dielectric 
function. Moreover, the method to obtain the parameter 
£ microscopically by considering the difference between 
the hypothetical bulk and the real equation of motion for 
surface degrees of freedom, can be generalized to other 
systems. In particular for optical phonons (polaritons) in 
insulators and photonic crystals (5^,^21 some of the 
above deviations from the conventional Fresnel theory 
can be expected. 
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tronic Properties-MaNEP" . 



APPENDIX A: EIGENMODES FOR 
ALTERNATING JUNCTIONS 

We introduce the unit cell, which contains two differ- 
ent junctions and describe the system by the parameters 
Cmi, dmi, Pmi (cf. Eqs. whcrc m denotes the unit cell 
and I — \,2 labels the junctions in the unit cell. The 
equations inside the crystal are analogous to Eqs. (^9|)- 
( pT| ) , where the quasiparticle dissipation is taken into ac- 
count by — ujf /uj'^Q I and ujf /lo'^ = 1 — Amaiuj / Ccoi^cO f 



^2?7 ^ - d,n2r] - c,n,iV + dmirj ^ = 0, 



(Al) 



2(Cmi77"^ - Cm-l^Vl) + (a - l)(P„ii - P,„-l,2) + 

i{P/b){cmi'>r' ~ dmiv) = 0, 

2(Cm2?7~^ - CrnlV) + (« - 1)(-Pm2 

i(/3/6)(c,„2?7~^ - d,n2r]) = 0, 

Pmliwi - a) + aiP,n2 + P7n-1,2 " 

(c,„i + rfmi)(l - 2af3){sinb/b), 

Pm2{w2 - a) + a{P,n+l,l + P„il - 2P„2) 

(c,„2 + rfm2)(l - 2a/3)(sin6/6). 



P„ 



IP 



(A2) 



(A3) 



Using the Fourier transformation with respect to the 
discrete index m we obtain the dispersion relation in the 
limit h <C <7,/3^/^; which is appropriate for oblique inci- 
dence in Sect. 

{y^ + 2(3- 0^)V + 4a(a - 1)(1 - i/^)/? +{a- ifv"^ + 

+/3)(a- l)(wi + W2 - 2a-4a) = 0, (A4) 
V^{wi~a- 2a)(w2 - a - 2a) - Aa^il - v^). 

In the opposite limit q <^h used for parallel incidence 
in Sect. ^ we get the dispersion 



Det 



We 



a-l 



0, 



2 + (3' 

hi = [l - t^o.i(l + 2a)/wf] , £12 = £21 = 2a, 



where W, i and A are matrices, Wu = wi, Wu = M^2i = 
0, and A;; = 1 + (3, A12 = A21 = 1. 

Using Eqs. (A4) for b <C q,P^^'^ and taking into ac- 
count the lattice constant 2s in qp — 2skz,p, the refraction 
indices Up = ckz,p/oJ — cvpj su: of the bulk eigenmodes 
for fj; = are determined by (zq = a — 1 — 2a) 



vl,2{^) = {-P±VP^ - 16a2/3g)/8a^ (A6) 
P{w) = w^S - wil + 6){l + 2a) + l + 4a(l - /?zo), 
Q{w) ^ 2w^6 - w{l + 5){1 + a + 4a) + 2a(l + 4a). 

Thereby and in the following the effect of dissipation can 
be included by replacing w ^ wi and wS — > W2 and 



we will restrict the discussion to the case a^ > (a — 
1)/3(1 - <5)/4 and na,/ 10,0,1^,0 < a[(a - l)/3]^/V[(l " 
(5)(1 + 2a)] when an extremal point We exists in the 
lower band, provided that 1 — (5 is of order unity. 

At a = 1 we get Q{w) = 2P{w) + 0{a^l3) and 
v1 2 is small near the zeros wiow,up of P(w), which are 
given by Eq. ([TSf). The reflection coefficient is deter- 
mined predominantly by small i'i2(w'), as in the case 
of identical junctions. Therefore, in the following we 
will analyze the behavior of Vi 2(10) by expanding around 



Wlow 

w — 



for the lower (upper) band. With u\o 



u'low.up we obtain P{w) = +\u\o 



and Q{w) 



-[±2A + (a - 1)(1 -I- 5)]wiow,up ± (a - l)iiow,up, where 
we denote A = [(1 - 5f{l + 4a) + Aa'^{l + SfY'"^ and 
iiow,up = ±[2 -I- 8a - wiow,up(l + 5)] > (upper/lower 
sign for lower/upper band). From this the band edges 
w[^2 up ^^"^ special frequencies We.i of the bands can be 
obtained in the limit ai = 0, cf. Fig. Ill] and [ij. 

In the lower band positive real solutions for exist 



for ai — aX We = w\ow 
where u\o^ = [16a^/3(a — 
point We we get — V2 = 

„(+) 



< w < w 

11/2 



(+) 
low 



= 1 2a, 



l)iiow]^"/A. At the extremal 
[/5(a-l)Liow]^/V2a, while the 

upper edge w\^J is determined by the condition = Ihy 
noting that P(l -|- 2a) = — 4a^. The value j^l is positive 



at w < wq 



+ ia- l)Liow/(2A + (a - 1)(1 + 6)) 



and approaches —2/3 for w > wq till the second band is 
reached. In the following we consider the case wl^2 > ""^o 
and hence the upper edge of the lower band is wl^2 ■ I'^i 
the range [ii'e,ii'o] two propagating modes with normal 
and anomalous dispersion exist, while for w G [wq, w^^J], 
vi is propagating and 1^2 decaying, which is very similar to 
a system with identical layers. Also note that the width 
Wo — We of the resonance in Tiw) is not proportional to 
the c-axis coupling a, but is mainly given by the angle of 
incidence. 

In contrast to this, the behaviour of and j^l in the 
upper band is quite different because the dispersion here 
is anomalous at any frequency. In this range is neg- 
(A5) ative and > 2(3 and the band edges are determined 
by the conditions i/J = or i/J = 1 respectively. The 



value v( is positive inside the band wip'' < w < w)^' , 
where 



„(+) 



„(-) 



= (1 + 2a)/ 5 and Wup'' = Wup 



Mup where 



^up 



= -{a - l)Lup/(-2A +{a- 1)(1 + 5)) > 0. At 



W = Wj 



= w^p we obtain — i/, = 



= = 



the point 

[P{a— l)Lup]^/^/2a, which corresponds to the frequency 
of maximal transmission according to Eq. ( po|) . 

Similarly to Eqs. (|52|)-(|5^) we make an ansatz for the 
bulk eigenvectors 



Cm/ 



= X! 7pC/(9p)cxp(igpm), 



p=l,2 



dmi= ^ 'ypdi{qp)exp{iqpm), 



(A7) 
(AS) 
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Pmi= ^ 'jpPi{qp)exp{iqpm). 

p=l,2 



(A9) 



Using Eqs. (|Al|)-( A3) we obtain the coefficients (P 
(Pi(g),P2(g))) 



ci = 1, C2(g) 

di = 



k(a- l)(l + e'"?) 



V+{a- 1){W2 - a - 2a)' 
l-[l + (772-,y-2)c2(q)]e-'« 



7/2 — Tj 



£ = M (c + rf) . 



M 



1 / W2 - a - 2a -a(l + e-''?) 
-a(l + e*«) wi-a-2a 



(AlO) 

(All) 

(A12) 
(A13) 
(A14) 



Here 02(9) is given in leading order in b and /3^/^. In 
order 0{iP) we get c?i = ci = 1 and = c^iq)- 

To determine 71/72 we use the microscopic boundary 
condition for the surface junction (analogous to Eq. (pO|)). 
Near We^i we get in leading order in (3^/"^ ~ ^a/n^ 

Poiiwi - a) + a{Po2 - [1 + a + (a - l)eeo]Poi} = 
(coi+doi)[l + a(ecO + l)]- (A15) 

Again we simplify this equation by subtracting the (hy- 
pothetical) bulk equation for Pqi, which follows from 
Eq. (A3) with Pm=-i.2 given by Eq. ([A9|), and the real 
surface Eq. ( |A15| ) for Pqi: 



Pm=-i,2 = ^7pP2(qp)exp(-i9p) = 0. 



(A16) 



Note that this ABC has only been derived in leading 
order in Ca/^-p a-nd near the resonance frequencies We^i- 



APPENDIX B: REFLECTIVITY IN PARALLEL 
INCIDENCE 

For arbitrary /3o and a the solutions of the dispersion 
Eq. ([A5|) in the case g <C & are given by: 



c^^xp _ {w6{w - Wpoio) + «<S'i)(l + /3o) 



(3o[w^S - 2awil + 6) + iS2] 
K 



1± 1 



(Bl) 

1/2 



{w5{w-Wp,y,)+iSiY{l + (3QY 

K = 2(3q{w'^ 5 - 2aw{l + 5) + iS2) (B2) 

{5{w - w\ow){w - w^p) + iS){l + Po/2) 
S ^ w^^^[{2a + l)6w{ai + ^2) - (1 + 4a){ai + ^26)], 

51 =u;3/2^(2a + l/2)(CTi+CT2), 

52 = 2aw3/2(5((Ti +CT2). 

Away from the pole lUpoic in Ec we obtain in leading order 
in Po for arbitrary a 



ai 



c _ 5{w- -wio„)(w - Wup) + iS 

Ld'^EcO W(5(W - Wpolo) + iS*! - /3oCo ' 

{w5{w - Wpoic) + iSi)[l + Po) 



a2 Lu^e^o {PQ/2){w^S~2aw{l + 6)+iS2)' 



(B3) 



(B4) 



In the ansatz Eq. (|77| ) for the field By outside the crys- 
tal the continuity equation at x = gives for By{kz) the 
following expression 



ckx 



dzG{z) exp{—ikzz), 



(B5) 



where for 2ms < z < (2m -I- l)s we obtain {gp = 

1/2 l/2lN 



G{z) 



c^^'' exp(i(7pz) -I- df exp{~igpz), (B6) 



j(p) 



p=l,2 



and for (2m + l)s < z < (2m -f 2)s analogously with 
c(f) _>4p) and ^''^ ^4"^. 

We derive — ±gp + {tt/s)], where j is an inte- 
ger. For nonzero j we obtain > tt/s and hence k^p 
is imaginary with large \kr^p\. For j = we obtain 
{B{gp)) = {B{—gp)) = by averaging over the two junc- 
tions in the unit cell, as c~f' + Sf' = d^f + djf = with 
accuracy 6//3o ~ Xab/Xc < 1 from Eqs. (|Ai|)-(|A^). As 
a result, in Eq. ( fn] ) the terms with amplitudes By{kz) 
may be dropped. Then the amplitude of the refiected 
wave By°^ is determined by Eq. (|^), where the averaged 
magnetic and electric fields at the boundary a; = are 



E ap[pi'^ + p^'^ 



(B7) 
(B8) 



p=l,2 



These equations lead to the reflection coefficient, Eq. (|^), 
where 



fli -I- a2Z 



Z = 



P'l 



(1) 



P 



(1) 



P 



(2) 



P 



(2) 



(B9) 



with rixp = ckxp/oJ. As 02 <C ai for w 7^ Wpoio the con- 
ventional Fresnel expression is valid everywhere except 
at the upper edge Wpoio of the lower band. 

To determine Z , we use the additional boundary con- 
dition P^ = of the Pekar type in the form 



Ex{x = 0) = ^ 



1 



1^2 ^xp^p 



1/2 



(v) (v) 

From this and the relation cth = dm the condition 
^pC^i"* /[kxpO'p] = follows. To express it in terms of 



P 



(p) 



P!f^ we derive from Eqs. (^-(^) 
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c(^) = -(l/4)(ap-l)p(^)(l-/p), 



fp = 



(p) 



3(P) 



(wi-l- 2a){2ap + /jp) ~ (cp - l)(ap + /jp) 
{ap - l){ap + Po) - 2a{2ap + Po) 



Finally we obtain 

(l-/i)(l + /2)(ai-l)fc,2a2 



Z = 



(l-/2)(l + /i)(a2-l)A;,iai' 
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